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ABSTRACT 


A  computer  model  of  an  aviation  gas  turbine  test 
facility  was  developed  and  analyzed  to  compare  computer 
generated  gas  pressures,  temperatures,  and  velocities  to 
those  of  an  actual,  operational  test  facility.  The 
computer  model  simulated  the  Royal  Danish  Air  Force  Test 
Facility  at  the  Skrydstrup  Air  Force  Base,  Denmark.  The 
Parabolic,  Hyperbolic,  or  Elliptic  Numerical  Integration 
Code  Series  (PHOENICS)  was  used.  The  Computer  Model  output 
compared  well  with  actual  data  obtained  from  the  Denmark 
test  facility.  Tabular  Data  and  Figures  are  presented. 
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I.  INTRODUCTION 


A.  BACKGROUND 

The  United  States  Navy  desires  to  provide  standardized 
test  facilities  for  the  testing  of  its  aviation  gas  turbine  ' 

jet  engines  (Ref  1) . 

Today's  high  performance  and  high  technology  jet 

I 

engines  require  accurate,  consistent,  and  controlled 
testing  to  assess  performance  throughout  the  jet  engine's 
entire  operational  range  Modern  day  test  facilities  must 

( 

incorporate  state  of  the  art  equipment  and  instrumentation 
such  that  jet  engines  may  be  tested  in  or  out  of  the 
aircraft,  around  the  clock,  and  year  round. 

I 

B.  NATURE  OF  THE  PROBLEM 

Presently,  75  Navy  and  Marine  Corp  activities  worldwide 
operate  81  Jet  Engine  Test  Facilities  (Ref  1) .  Of  these, 

60  are  fully  operational  and  51  are  over  25  years  old. 

Immediate  problem  areas  which  hinder  operational  needs 
include;  structural  degradation,  obsolescence,  noise,  and 
air  pollution. 

1 .  Structural  Degradation 

High  jet  engine  exhaust  gas  temperatures  and  water 
impingement  on  diffusers  and  interior  concrete  walls  have 
caused  cracking  and  spalling  of  test  facility  interiors. 
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Concrete  reinforcing  rods  exposed  to  a  hot,  moist 
environment  will  corrode  and  thereby  reduce  structural 
integrity.  Noise  reduction  baffles  and  exhaust  gas  turning 
vanes  are  being  damaged  by  flying  concrete  debris. 
Replacement  of  baffling  and  turning  vanes  is  required  as 
often  as  every  three  years.  The  maintenance  requirements 
navy  wide  has  significantly  increased  the  workload  and 
operational  costs  of  the  Navy  Public  Works  Department. 

2 .  Obsolescence 

Several  of  the  test  facilities  require  major 
modernization  or  even  complete  replacement  due  to  current 
aircraft  designs  and  test  requirements.  Examples  of 
modernization  include  improved  engine  handling  systems  and 
engine  thrust  stands.  Increased  access  to  engine 
components  and  test  points  greatly  improves  testing 
efficiency.  Additionally,  with  today's  high  technology 
engines,  state  of  the  art  data  acquisition  systems  are 
required . 

3 .  Noise 

Most  of  the  older  test  facilities  use  water  for 
exhaust  gas  cooling.  The  water  and  the  high  exhaust  gas 
temperature  deteriorates  the  acoutic  baffling  silencing 
capabilities.  The  increased  .  noise  level  degrades  the  work 
environment.  Worker  efficiency  is  reduced  and  the  high 
test  facility  interior  sound  levels  are  being  suspected  of 
contributing  to  engine  deterioration.  State  of  the  art 


engines  using  thinner,  lighter  materials  are  more  sensitive 
to  acoustic  fatigue  (Ref  1) . 

4 .  Air  Pollution 

Today's  high  power  and  high  mass  flow  engines 
constitute  a  major  problem  to  the  environment.  Visible 
particulate  matter  is  emitted  in  excess  of  allowable 
environmental  pollution  standards.  By  Executive  Order 
Number  11282,  May  26,  1966  (Ref  2),  Jet  Engine  Test 
Facilities  must  meet  local  laws  concerning  the  environment. 
Judicial  action  by  the  state  of  California  against  Jet 
Engine  Test  Facilities  within  California  (Ref  3)  has 
identified  a  need  to  address  pollution  control  measures 
during  the  design,  construction,  and  operation  of  Jet 
Engine  Test  Facilities. 

C.  TYPES  OF  TEST  FACILITIES 

The  Navy  operates  five  types  of  Jet  Engine  Test 
Facilities  (Ref.  1).  Each  facility  utilizes  some  form  of 
hold  down  device  for  the  engine  or  the  aircraft  and 
instrumentation  to  monitor  the  performance  of  the  engine. 
Each  facility  must  be  designed  such  that  several  engine 
types  may  be  tested  (e.g.,  turboprop,  turbofan,  turboshaft, 
turbojet) .  The  types  of  facilities  are; 

1 .  Run-up  Pad 

The  simplest  test  facility  is  the  run-up  (power 


check)  pad.  Normally,  the  engine  is  tested  while  installed 
in  the  aircraft.  Therefore,  the  aircraft  acts  as  the  hold- 


down  device.  The  run-up  pad  is  usually  located  in  remote 
locations  due  to  high  noise  levels  encountered  during 
testing. 

2 .  Sound  Suppressors 

These  are  usually  equipment  packages  which  are  used 
to  reduce  inlet  and  exhaust  noise  from  a  tied  down 
aircraft.  The  sound  suppressors  may  be  used  to  augment  or 
to  substitute  for  the  run-up  pads. 

3 .  Open  Test  Stand 

This  is  an  outdoor  test  facility  which  supports  and 
constrains  the  engine  during  test.  Due  to  noise  levels, 
the  open  test  stand  must  be  located  in  a  remote  area. 
Portable  sensors  and  instrumentation  are  utilized  to 
monitor  engine  performance. 

4 .  Hush  House 

This  is  a  total  aircraft  enclosure.  It  permits 
round  the  clock  and  year  round  testing.  Normally  on  twin 
engine  aircraft,  only  one  engine  may  be  tested  in  full 
after  burner.  The  other  engine  is  held  at  idle. 

5 .  Fully  Enclosed  Test  Cell 

This  test  facility  is  normally  capable  of  testing 
several  types  and  sizes  of  jet  engines.  The  engines  may  be 
in  or  out  of  the  aircraft.  The  facility  is  normally 
equipped  with  all  the  necessary  instrumentation  and  data 
gathering  equipment  required  to  perform  routine  or  complex 
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testing  after  jet  engine  rework  or  overhaul.  Testing  may 
be  performed  round  the  clock  and  year  round. 

The  type  of  test  facility  to  utilize  at  an  activity  is 
dependent  upon  several  factors.  Availability  of  military 
construction  funds,  required  control  of  noise  and  emissions 
(environmental  concern) ,  amount  of  maintenance  testing,  and 
military  and  civilian  manpower  requirements  are  all  key 
factors. 

D.  RESEARCH  GOAL 

To  provide  quality  jet  engine  testing,  to  comply  with 
local  environmental  standards,  and  to  ensure  fleet 
readiness,  the  Navy  embarked  on  a  comprehensive  program  in 
1984  to  upgrade  and  to  modernize  as  necessary  all  jet 
engine  test  facilities  under  its  cognizance.  To  this  end, 
it  was  necessary  to  review  existing  jet  engine  test 
procedures  and  to  examine  jet  engine  test  facilities  to 
determine  which  facilities  met  desired  criteria. 
Additionally,  the  Navy  decided  to  explore  the  feasibility 
of  the  computer  modeling  of  aviation  gas  turbine  jet  engine 
test  facilities.  It  was  felt  that  if  a  computer  model 
could  accurately  represent  an  actual  operational  test 
facility  then  a  computer  model  could  also  be  used  to 
predict  the  operation  of  future  proposed  test  facilities. 
The  purpose  of  this  research  was  to  computer  model  the 
Royal  Danish  Test  Facility  and  to  compare  computer 


generated  temperatures,  pressures,  and  velocities  with 
actual  data  obtained  in  the  field  and  to  investigate 
techniques  necessary  to  make  computer  simulations  of  flow 
and  temperature  fields  useful  for  predictive  purposes. 


II.  TEST  FACILITY  OPERATION 


The  basic  geometry  of  the  Denmark  test  facility  is 
shown  in  Figure  C.l.  For  simplicity  all  storage  areas, 
engine  preparation  rooms,  and  control  modules  are  not 
included.  The  test  facility  is  basically  comprised  of  the 
primary  air  inlet,  a  test  bay,  an  exhaust  duct  (also  called 
an  augmentor  tube) ,  secondary  air  inlets,  and  an  exhaust 
stack.  Turning  vanes,  which  are  not  shown,  are  located  in 
the  exhaust  duct  to  exhaust  stack  transition  region.  This 
feature  better  utilizes  the  total  exhaust  stack  cross 
section  with  an  attendant  reduction  in  exhaust  gas  noise 
(Ref  4) .  Baffles  (not  shown)  are  located  at  three  places; 
At  the  top  of  the  exhaust  stack,  within  the  secondary  air 
inlets,  and  within  the  primary  air  inlet  (test  bay  doors) . 
The  baffles  also  provide  noise  reduction  (Ref  4). 

During  a  jet  engine  test,  the  engine  is  aligned  with 
the  centerline  of  the  exhaust  duct.  This  can  be 
accomplished  either  with  the  engine  in  the  aircraft  or  with 
the  engine  supported  in  a  test  stand.  With  proper 
alignment,  a  near  uniform  flow  profile  develops  (Ref  3) 
inside  the  augmentor  tube.  Primary  air  is  drawn  in  through 
the  primary  air  inlet  (test  bay  door/baffle  configuration) . 
For  a  F-lOO  Engine,  at  full  throttle  and  after  burner  on, 
the  primary  air  mass  flow  rate  is  260  kg/sec  of  which 
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100  kg/sec  is  the  mass  flow  rate  through  the  engine 
(Ref.  4) .  As  the  exhaust  gas  enters  the  augmentor  tube, 
air  from  within  the  test  bay  is  entrained  into  the 
augmentor  tube.  The  exhaust  gas  moving  through  the 
augmentor  tube  also  entrains  air  coming  through  the 
secondary  air  inlets.  The  secondary  air  mass  flow  rate  is 
340  kg/sec  (Ref.  4) .  Both  the  entrained  primary  and 
secondary  air  play  an  important  role  in  the  test  facility 
operation.  Hot  and  noisy  jet  engine  gases  must  be  treated 
prior  to  leaving  the  test  facility  in  order  to  comply  with 
environmental  standards.  To  treat  jet  engine  exhaust,  its 
temperature,  noise  level,  and  particulate  density  must  be 
lowered. 

The  ratio  of  total  "cooling”  air  (primary  and 
secondary)  mass  flow  rate  to  the  engine  mass  flow  rate  is 
called  the  Augmentation  Ratio.  The  Augmentation  Ratio  is  a 
key  to  successful  test  facility  operation.  With 
appropriate  primary  air  and  secondary  air  mass  flow  rates, 
the  exhaust  gas  is  sufficiently  cooled  which  helps  prolong 
duct/stack  life  and  the  exhaust  gas  is  sufficiently  diluted 
such  that  proper  pollution  levels  are  obtained  (Ref.  2) . 

The  Augmentation  Ratio  of  the  Denmark  Test  Facility  is  5:1. 

With  a  computer  model  that  can  accurately  predict  test 
facility  operation,  design  personnel  will  have  the  ability 
to  easily  change  model  parameters  such  as  engine  location 
or  structural  dimensions  and  to  observe  their  effect  on 
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test  facility  operation.  For  example  by  changing  the 
secondary  air  inlet  configuration,  the  augmentation  ratio 
is  changed.  This  in  turn  has  an  effect  over  exhaust  gas 
cooling.  Additionally,  design  engineers  may  wish  to 
experiment  with  various  exhaust  stack  heights  or  exhaust 
duct  openings.  A  mere  changing  of  a  few  lines  of  computer 
model  source  code  accomplishes  this.  In  all,  a  computer 
model  could  save  significant  time  and  money  in  the  design 
of  future  jet  engine  test  facilities. 


III.  DESCRIPTION  OF  PHOENICS 

A.  WHAT  IS  PHOENICS 

PHOENICS  is  a  general-purpose,  computer-code  system  for 
simulating  fluid  flow,  heat  transfer,  and  chemical  reaction 
processes  in  engineering  equipment  or  the  environment. 
PHOENICS  is  an  acronym  for  Parabolic,  Hyperbolic,  or 
Elliptic  Numerical  Integration  Code  Series.  Typical 
applications  of  PHOENICS  include:  Thermal  hydraulics  of 
nuclear  reactors  and  other  power  plants;  external 
hydrodynamics  and  aerodynamics  of  cars,  ships,  buildings, 
etc;  air  movement  within  an  enclosed  space,  both  under 
normal  conditions  and  in  the  presence  of  a  fire;  flows  in 
pumps  and  turbomachinery.  PHOENICS  was  developed  by  D. 
Brian  Spalding  of  CHAM  Ltd  and  of  Imperial  College  of 
London,  England.  Information  regarding  the  history  and 
development  of  PHOENICS  is  contained  in  the  user's  guide 
(Ref  5) . 

B.  HOW  PHOENICS  WORKS 

The  PHOENICS  computer  code  consists  of  two  separate 
computer  programs;  EARTH  and  Ql.  The  EARTH  program 
resides  in  the  computer  mainframe  as  an  object  code  and 
cannot  be  accessed  by  the  user.  The  user  couples  his 
problem  (defined  by  the  user)  to  the  EARTH  code.  The  EARTH 


10 


program  contains  the  main  flow  simulating  software  and 
various  coding  sequences  which  simulate  the  Laws  of  Physics 
as  applied  to  the  conservation  principles:  energy,  mass, 
momentum.  The  Q1  program  defines  the  problem  and  differs 
from  one  flow  simulation  program  to  the  next.  The  Q1 
program  represents  the  user's  decision  as  to  what  problem 
is  to  be  simulated.  The  Q1  program  activates  the  EARTH 
program  which  solves  the  problem. 

The  PHOENICS  code  is  capable  of  describing  flow 
simulations  in  terms  of  distributions  in  time  and  space. 

The  distributions  involve  assigning  numerical  values  for 
all  dependent  variables  to  an  orderly  array  of  cells 
located  in  the  physical  domain  of  the  problem. 

Temperature,  pressure,  concentration,  KE  (Turbulence 
Kinetic  Energy) ,  and  EP  (Kinetic  Energy  Dissipation  Rate) 
values  are  located  at  a  cell's  center  while  velocity  values 
are  located  at  the  positive  faces  of  the  cells.  A  cell 
configuration  in  the  X-Y  plane  is  illustrated  in  Figure 
C.2.  Temperatures,  pressures,  and  concentrations  would  be 
located  at  N,  S,  E,  W,  and  G  for  example  while  velocities 
would  be  at  N',  S',  E',  and  W' .  A  similar  procedure  is 
utilized  for  the  Z  (low-to-high)  direction.  This 
convention  is  called  the  "STAGGERED  GRID"  arrangement.  If 
the  flow  simulation  is  time  dependent,  the  distributions 
are  obtained  for  each  variable  at  a  succession  of  times 
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instants.  The  Denmark  test  facility  model  flow  simulation 
is  for  steady  state  conditions. 

In  general,  a  PHOENICS  flow  simulation  is  four 
dimensional:  1  time  and  three  space  coordinates.  For  ease 

of  reference  to  the  user,  PHOENICS  measures  time  in  the 
past-to-future  direction;  X  coordinates  in  the  west-to-east 
direction,  Y  coordinates  in  the  south-to-north  direction, 
and  Z  coordinates  in  the  low-to-high  direction.  This 
concept  is  illustrated  in  Figure  C.3.  Although  Figure  C.3 
represents  a  Cartesian  configuration,  individual  cells  may 
vary  in  shape  and  size  as  determined  by  a  user  defined 
code.  This  option,  called  BFC,  was  not  utilized  in  this 
investigation.  The  Denmark  test  facility  model  is  a 
Cartesian  configuration  of  cells  which  vary  in  size. 

PHOENICS  calculates  the  dependant  variables  over  the 
complete  set  of  cells,  which  defines  the  operational  (e.g., 
system,  physical)  domain,  with  linear  algebraic  equations 
that  represent  the  conservation  laws.  Here  the  partial 
differential  equations  are  those  governing  the  conservation 
of  energy,  momentum,  and  mass  and  the  linear  algebraic 
equations  are  obtained  by  the  finite  volume  method.  The 
finite  volume  method  is  a  method  having  some  features  from 
the  finite  differencing  method  and  the  finite  element 
method.  The  conservation  partial  differential  equations 
will  not  be  discussed  here;  except  to  say  that  the  original 
non  linear,  coupled  conservation  partial  differential 


equations  are  linearized.  Interested  readers  are  directed 
to  the  PHOENICS  User's  Manual  (Ref.  5)  for  detailed 
information. 

PHOENICS  utilizes  an  iterative  process  involving  a 
multi-stage  sequence  of  dependent  variable  value 
adjustments.  The  process  is  repeated  several  times  until 
the  left  and  right  imbalance  (called  a  residual)  of  each 
equation  is  reduced  to  a  negligible  amount.  A  variety  of 
iterative  procedures  are  available  to  the  user.  The 
"SLAB-SWEEP"  iterative  process  feature  of  PHOENICS  for 
equation  imbalance  reduction  was  chosen  for  the  Denmark 
test  facility  model.  Slabs  are  arrays  of  cells  having  the 
same  Z  (low-to-high)  coordinate.  PHOENICS  conducts  many 
mathematical  operations  over  a  slab  in  an  attempt  to  reduce 
dependent  variable  imbalances  among  cells  within  the  slab 
prior  to  moving  on  to  the  next  slab.  The  adjustment 
operations  are  referred  to  as  "SLABWISE  SOLUTIONS".  The 
operations  begin  at  the  lowest  numbered  slab  and  continue 
sequentially  to  the  highest  numbered  slab.  A  sweep  is  a 
completed  set  of  slabwise  operations.  The  equations  for 
dependent  variable  values  at  one  slab  are  coupled  to  values 
of  other  dependent  variables  at  neighbor  slabs.  For  this 
reason  numerous  sweeps  are  normally  required  until  all 
dependent  variable  imbalances  are  negligible.  This  is  a 
necessary  but  not  a  sufficient  condition  for  problem 


convergence . 


C.  PHOENICS  INPUT  LANGUAGE 


The  source  code  used  in  a  PHOENICS  Q1  program  is 
arranged  in  a  24  group  data  structure.  Figure  C.4  lists 
the  group  categories  of  the  data  structure.  The  developers 
of  the  PHOENICS  code  suggest  this  ordering  of  the  groups 
normally  conforms  to  a  well  thought  out  flow  simulation 
problem  description.  However,  any  group  ordering  the  user 
desires  will  be  acceptable  to  PHOENICS.  It  is  the  PHOENICS 
defined  variables  and  keywords  plus  the  user  defined 
variables  that  actually  trigger  desired  PHOENICS  responses 
as  opposed  to  a  particular  group  ordering.  A  detailed 
description  of  group  functions  is  explained  in  the  PHOENICS 
User's  Manual  (Ref.  5) . 

D.  MODEL  GEOMETRY 

PHOENICS  utilizes  three  distinct  procedures  for 
defining  a  flow  simulation  operational  domain  (grid) : 
Cylindrical-polar,  Body-fitted,  and  Cartesian  coordinate 
systems . 

1.  Cvlindrical-Polar  Coordinate  System 

The  cylindrical-polar  grid  arrangement  consists  of 
planes  with  constant  Z  coordinates  that  are  perpendicular 
to  an  axis  of  rotation,  planes  with  a  constant  X  coordinate 
(usually  called  the  Theta  coordinate)  that  intersect  with 
the  axis  of  rotation  giving  the  X  coordinate  a 


representation  in  radians,  and  concentric  cylindrical 
surfaces  with  a  constant  Y  radial  coordinate. 


2 .  Body-Fitted  Coordinate  System 

This  arrangement  is  best  thought  of  as  starting 
with  an  orthogonal  group  of  cells  made  of  clay.  The  cells 
are  then  squeezed,  bent,  twisted,  and  stretched  into  the 
desired  shape  for  the  fluid  flow  simulation  operational 
domain.  Each  cell  although  distorted  in  shape  and  size 
still  remains  in  contact  with  its  original  neighbors. 

3 .  Cartesian  Coordinate  System 

The  Cartesian  grid  arrangement  is  composed  of  cells 
in  a  mutually  perpendicular  X,  Y  and  Z  coordinate  system. 

As  such  the  cell  faces  form  planes  that  intersect  at  right 
angles.  The  spacing  between  the  cell  faces  (planes)  is 
user  defined. 

Due  to  the  regularity  of  a  Cartesian  coordinate 
system,  the  information  required  to  specify  the  geometry  is 
rather  modest.  The  PHOENICS  yariables  NX,  NY,  and  NZ 
specify  the  number  of  cells  the  operational  domain  is  to 
haye  in  the  X,  Y,  and  Z  directions  respect iyely.  The 
PHOENICS  yariables  XULAST,  YVLAST,  and  ZWLAST  specify  the 
domain  lengths  in  the  X,  Y,  and  Z  directions.  The  "METHOD- 
OF-PAIRS"  is  a  PHOENICS  coding  sequence  to  specify  the 
number  of  cells  per  length  of  domain  in  the  X,  Y,  and  Z 
directions.  XFRAC,  YFRAC,  and  ZFRAC  are  PHOENICS  yariables 
for  arrays  used  in  the  method-of-pairs  coding  sequence. 


The  odd  indexed  array  elements,  for  example,  XFRAC(l), 
XFRAC(3),  XFRAC(5),  etc.,  specify  the  number  of  cells  in  a 
particular  length  of  operational  domain  while  the  even 
indexed  array  elements  XFRAC(2) ,  XFRAC(4),  and  XFRAC(6) 
specifies  the  length  of  each  cell.  A  negative  argument  for 
the  first  element  in  each  array  tells  PHOENICS  that  the 
method-of-pairs  is  being  utilized.  For  example,  ''XFRAC(l) 
=-2.0;XFRAC(2)=0.1''  tells  PHOENICS  the  method-of-pairs  is 
to  be  used  and  that  the  first  cell  group  in  the  X  direction 
is  to  have  two  cells  and  the  length  of  each  cell  is  0.1 
meter.  "XFRAC(9)=4.0;XFRAC(10)=0.15"  specifies  the  number 
of  cells  and  lengths  of  each  cell  in  the  "FIFTH"  X 
direction  cell  group. 

E.  POROSITY  AND  BLOCKAGE 

In  any  flow  simulation  problem  there  exists  walls, 
submerged  bodies,  or  some  form  of  obstruction  to  fluid 
flow.  For  example,  in  Hagen-Poiseuille  pipe  flow  the  fluid 
is  confined  to  the  interior  of  the  pipe.  Additionally,  due 
to  wall  friction  the  fluid's  velocity  at  the  wall  is  zero. 
In  a  PHOENICS  flow  simulation  problem,  obstructions  to 
fluid  flow  must  also  be  introduced.  In  a  jet  engine  test 
facility  model  one  must  be  concerned  with  obstructions  to 
flow  as  a  result  of  walls,  floors,  baffle  material,  turning 
vanes,  and  the  jet  engine  itself.  To  account  for  the 
physical  obstructions,  PHOENICS  makes  use  of  the  command 
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"CONPOR(R.R,  TYPE,XF,XL,YF,YL,ZF,ZL) The  real  number 


"R.R"  normally  has  a  value  from  0.0  to  1.0.  This  value 
tells  PHOENICS  how  much  obstruction  to  flow  is  to  be 
introduced  by  a  surface  or  volume.  For  maximum  obstruction 
(blocked  flow),  "R.R”  would  be  0.0.  The  “TYPE"  argument 
specifies  a  surface  or  volume.  For  example  a  surface  might 
be  the  inside  of  the  test  bay  ceiling.  A  volume  might  be 
the  geometric  body  of  the  jet  engine.  For  a  surface 
obstruction  the  "TYPE"  Arguments  "WEST,  EAST,  SOUTH,  NORTH, 
LOW,  OR  HIGH"  (the  faces  of  an  individual  cell)  are  used. 
For  a  volume  obstruction  the  type  "CELL"  is  used.  The 
arguments  "XF,  XL,  YF,  YL,  ZF,  and  ZL"  identify  the  first 
and  last  cells  of  the  obstruction  domain  in  the  X,  Y,  and  Z 
directions  respectively.  For  example  "CONPOR(0.0,  HIGH, 
1,3, 1,1, 3, 3)"  tells  PHOENICS  to  introduce  complete 
obstruction  to  flow  through  the  "HIGH"  (Z  direction) 
surfaces  of  the  group  of  cells  starting  at  IX=1,  IY=1,  and 
IZ=3,  and  ending  at  cells  IX=3,  IY=1,  and  IZ=3 
respectively.  Each  cell  has  an  IX,  lY,  and  IZ  identity. 
Figure  C.5  illustrates  "HIGH"  surface  obstruction  to  flow. 
The  cross-hatched  surface  blocks  Z  direction  flow.  "CONPOR 
( 0 . 0 , CELL, 1 , 3 , 1 , 1 , 1 , 3 ) "  would  introduce  maximum  obstruction 
to  flow  through  the  surfaces  of  the  volume  defined  by  these 
coordinates.  In  Figure  C.5,  the  "CELL"  volume  outlined 
with  bold  lines  prevents  fluid  flow  perpendicular  to  any  of 


its  surfaces. 


PHOENICS  also  has  the  ability  to  introduce  wall  or 
surface  "FRICTION"  which  brings  about  zero  magnitude 
velocity  components  on  a  cell's  surface  and  calculates  the 
friction  coefficient.  For  the  above  "HIGH"  example,  the 
magnitude  of  flow  in  the  Z  direction  was  set  to  zero  on  the 
surfaces  of  several  cells.  To  also  set  the  X  and  Y 
velocity  components  to  zero  on  the  same  "HIGH"  surfaces, 
the  CONPOR  command  requires  a  slight  modification. 

Utilizing  a  "NEGATIVE"  argument  with  the  CONPOR  command 
accomplishes  the  desired  action.  For  example 
"CONPOR(0.0,HIGH,l,3,l,l,3,-3)"  causes  the  X,  Y,  and  Z 
velocity  components  on  the  "HIGH"  surfaces  of  that 
particular  group  of  cells  to  be  zeroed.  Through  the  use  of 
the  negative  argument  convention,  velocity  components  on 
the  external  surfaces  of  a  volume  (cell)  can  be  zeroed,  and 
friction  introduced  into  the  model . 

F.  SOURCE  AND  BOUNDARY  CONDITIONS 

In  any  PHOENICS  flow  simulation  problem  description, 
the  user's  Q1  program  must  introduce  source  and  boundary 
conditions  to  the  EARTH  program.  Through  the  use  of 
PHOENICS  keywords,  source  and  boundary  conditions  may  be 
identified  as  energy,  momentum,  or  mass  flow.  Energy 
sources  are  identified  with  the  enthalpy  variable  "H", 
momentum  by  the  velocity  variables  "U,  V,  or  W" ,  and  mass 
flow  utilizes  the  pressure  variable  "P".  By  using  the 


notation  ''HI”  a  fluid's  first  phase  enthalpy  would  be 
identified  while  ''H2''  would  refer  to  second  phase  enthalpy 
for  a  two  phase  flow  simulation.  Information  about  inflow 
and  outflow  of  any  variable  in  the  system  domain  is 
conveyed  to  the  EARTH  Program  by  way  of  the  PHOENICS 
combination  "PATCH-COVAL”  statements.  The  PATCH  command  is 
used  to  specify  the  location  at  which  a  source  or  boundary 
condition  is  to  be  applied  and  to  also  specify  the  "TYPE” 
factor  or  what  the  PATCH  command  is  to  do  to  the  system 
domain.  The  COVAL  command  is  used  to  specify  the 
coefficient  and  value  which  combine  to  determine  the 
magnitude  of  the  dependent  variable  over  a  PATCH.  For 
example,  "PATCH(JETOUT,LOW,l,l,5,6,8,8,l,l) ; 

COVAL(JETOUT, PI, FIXFLU, 100.0) ; 

COVAL (JETOUT , W1 , ONLYMS ,500.0); COVAL (JETOUT , HI , ONLYMS , 
200000.0)"  specifies  a  100  kg/sec,  steady-state,  fixed-flux 
mass  source.  The  mass  source,  named  JETOUT,  begins  at  the 
low  faces  of  the  cells  IX=1,  IY=5,  and  IZ=8  and  ends  at 
IX=1,  IY=6,  and  IZ=8.  The  mass  source  has  a  Z  direction 
velocity  of  500  m/sec  and  an  enthalpy  of  200,000  J/kg.  The 
COVAL  coefficient  "ONLYMS"  specifies  that  only  convection 
flow  effects  are  influential  in  the  transport  of  enthalpy 
from  the  cell  faces  to  the  cell  centers.  That  is, 
diffusive  effects  are  zero.  Appendix  2  of  the  PHOENICS 
User's  Manual  (Ref.  5)  lists  the  various  "TYPES", 


"COEFFICIENTS”,  and  "VALUES"  that  may  be  used  with  the 
combination  "PATCH-COVAL"  statements. 

G .  INITIAL  CONDITIONS 

In  order  to  obtain  the  solution  of  a  steady  state 
problem,  PHOENICS  starts  the  iterative  solution  with  an 
initial  field  of  the  dependent  variables.  The  default 
initial  fields  for  all  dependent  variables  is  zero  over  the 
entire  physical  domain.  However,  convergence  to  the  true 
solution  can  be  speeded  up  by  the  user  specifying  initial 
field  conditions  via  "FIINIT"  (Field  Initialization) 
statements  over  the  entire  operational  domain  or 
"PATCH-INIT"  statements  over  subdomains  of  the  system.  For 
example  "FIINIT (HI) =AMOUNT"  initialize  a  first  phase 
enthalpy  field  to  a  user-defined  value  "AMOUNT". 

"PATCH (GASIN , INIVAL, XF, XL, YF, YL, ZF, ZL, TF , TL) ; INIT (GASIN , PI , 
COEFFICIENT,  VALUE)"  initializes  a  first  phase  pressure 
field  for  the  user-defined  PATCH  "GASIN".  XF,XL, YF, YL, ZF, 
and  ZL  define  the  geometric  region  of  the  PATCH  called 
"GASIN". 

H .  MAXIMUM/MINIMUM  VALUES 

At  the  beginning  of  a  flow  simulation  problem  iterative 
process,  a  poorly  defined  dependent  variable  field  value 
may  have  an  adverse  affect  on  problem  convergence.  At 
times,  dependent  variables  "STRAY"  outside  of  a  physically 
meaningful  range  and  cause  numerical  divergence.  To  help 
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minimize  divergence  of  the  solution  PHOENICS  utilizes  the 
"VARMIN"  and  VARMAX"  commands.  These  commands  are  used  to 
set  upper  and  lower  bounds  on  any  dependent  variables. 
"VARMIN” (VI) =SLOW; VARMAX (W1)=FAST"  sets  a  lower  limit  of 
"SLOW"  to  VI  and  an  upper  limit  of  "FAST"  to  Wl. 

I .  UNDER-RELAXATION  FACTORS 

A  feature  of  PHOENICS  which  greatly  assists  problem 
convergence  is  the  application  of  "UNDER-RELAXATION" 
factors  to  flow  simulation  dependent  variables.  This 
practice  is  most  useful  on  an  operational  domain  comprised 
of  ill-proportioned  cells  (very  large  cells  adjacent  to 
very  small  cells) .  During  program  execution,  EARTH'S 
dependent  variable  correction  equations  may  develop 
enormously  large  coefficients  between  adjacent  cells.  The 
coefficient  disparity  often  contributes  to  solution 
divergence.  With  selective  "UNDER-RELAXATION"  applied  to  a 
flow  simulation  program,  PHOENICS  introducers  "FALSE-TIME- 
STEP"  values  to  select  dependant  variables  thereby 
preventing  large  dependent  variable  changes  in  magnitude 
per  sweep.  The  "FALSE-TIME-STEP"  brings  equation 
coefficients  closer  together  "LESS  ABRUPTLY"  which  promotes 
problem  convergence.  "UNDER-RELAXATION"  for  pressure, 
velocity,  and  enthalpy  components  is  frequently  required 
for  steady  state  problems. 
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J.  PROGRAM  RUN  PROCEDURES 

The  Denmark  test  facility  model  was  analyzed  on  the 
Naval  Postgraduate  School  IBM-3033  mainframe  computer 
system.  Program  run  procedures  involve  interactive 
terminal  sessions  and  a  batch  submission.  The  necessary 
flow  simulation  Q1  program  and  Job  Control  Language  (JCL) 
programs  reside  in  the  User's  mainframe  storage  area.  All 
flow  simulation  programs  are  submitted  for  execution  by 
means  of  "EXEC”  (execute)  and  JCL  files.  The  EXEC  files 
are  activated  by  the  keywords  "COPYQl",  "RUNSAT”,  and 
"RUNEAR" . 

The  COPYQl  EXEC  file  queries  the  user  as  to  the  name 
and  file  type  of  the  flow  simulation  Q1  program  source 
code.  Afterwards  the  COPYQl  EXEC  file  stores  the  source 
code  and  appropriate  JCL  in  a  temporary  file.  At  this 
point,  the  user  must  ensure  the  above  transfer  of  data  was 
successful.  This  is  accomplished  through  the  user's 
"VIRTUAL  READER"  file.  An  inspection  of  this  file  will 
reveal  various  condition  codes  which  indicate  the  status  of 
data  transfer.  Figure  C.6  is  an  example  of  COPYQl  EXEC  and 
JCL  files. 

The  "RUNSAT"  EXEC  file  extracts  from  the  user's 
mainframe  storage  area  JCL  which  enables  a  Q1  program  to 
EARTH  program  communications  link.  This  JCL  is  also  stored 
in  a  temporary  file.  Again  a  virtual  reader  file 
inspection  must  be  made  to  ensure  that  a  successful  data 


transfer  was  achieved.  Figures  C.7  and  C.8  are  examples  of 
RUNSAT  EXEC  and  JCL  files. 

Finally  after  successful  COPYQl  and  RUNSAT  data 
transfers,  the  flow  simulation  problem  is  ready  for 
execution.  The  "RUNEAR”  EXEC  file  is  called.  Necessary 
mainframe  storage  is  allocated  and  the  PHOENICS  system  is 
activated.  Figure  C.9  is  an  example  of  RUNEAR  EXEC  and  JCL 
files.  Upon  program  completion,  PHOENICS  output  is  sent  to 
the  user's  virtual  reader  file.  After  the  PHOENICS  output 
is  analyzed,  the  user  may  discard  the  information  or 
utilize  a  PHOENICS  restart  which  allows  the  problem  to 
continue  from  where  it  last  ended.  Upon  completion  of  a 
PHOENICS  run,  output  is  stored  on  disk  file  ''DF09".  To 
access  this  file  and  to  continue  the  problem  the  PHOENICS 
code  "RESTRT(A,  B,  .  .  . . ,  Z)  •'  must  be  utilized  in  the  Q1 
program.  ”A,B,...,Z”  represent  the  dependent  variables 
which  are  to  continue  from  their  last  value  when  the 
program  ended.  If  the  user  desires  all  dependent  variables 
to  restart  from  their  last  values  then  "RESTRT (ALL) "  can  be 
used  in  the  Q1  program. 


IV.  MODELING  THE  DENMARK  TEST  FACILITY 

A.  MODEL  DESCRIPTION 

The  Denmark  test  facility  model  utilizes  a  Cartesian 
grid  arrangement.  Figure  C.IO  is  a  schematic  diagram  of 
the  test  facility  in  the  Y-Z  pane  (side  view).  Figure  C.ll 
is  the  X-Z  plane  (top  view)  representation.  The  PHOENICS 
code  allows  for  a  plane  of  symmetry.  For  the  Denmark  test 
facility  model  the  vertical  Y-Z  plane  is  one  of  symmetry. 
Therefore,  in  Figure  C.ll  only  half  of  the  Denmark  test 
facility  is  represented.  In  this  figure  the  jet  engine  is 
slightly  off-of-center  from  the  Z  axis  for  clarity  only. 

To  facilitate  the  Denmark  test  facility  model  geometry 
specifications,  the  system  domain  was  divided  into  user 
defined  zones.  Each  zone  comprises  a  particular  section  of 
the  Denmark  test  facility  model.  For  example  in  Figure 
C.12,  "YLBLl”  is  a  user  defined  variable  which  specified 
the  length  of  the  first  zone  in  the  Y  direction.  This 
particular  zone  represents  the  length  in  the  Y  direction 
from  the  ground  or  floor  of  the  Denmark  test  facility  up 
2.0  meters.  "ZLEXT2"  represents  the  length  in  the  Z 
direction  of  the  atmosphere  behind  the  test  facility,  and 
so  on.  Figure  C.12  is  a  schematic  diagram  of  the  system 
domain  zones  which  make  up  the  vertical  Y-Z  plane  while 
C.13  illustrates  the  zones  in  the  horizontal  X-Z  plane. 
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The  zone  concept  has  proven  to  be  very  flexible.  The 
number  of  cells  in  each  zone  may  be  specified  by  a  user 
defined  variable.  For  example,  *'XNCBL1''  specifies  the 
number  of  cells  in  the  X  direction  for  the  length  of 
operational  domain  ''XLBLl”.  *'ZNCEX2''  specifies  the  number 
of  cells  in  the  Z  direction  for  the  length  of  operational 
domain  ''ZLEXT2"  and  so  on-  Grid  refinement  for  a 
particular  zone  can  be  accomplished  without  any  difficulty 
by  changing  one  variable  in  the  Q1  code.  Figure  C.14  is 
the  method-of-pairs  coding  sequence  which  specifies  the 
operational  domain  zones  and  numbers  of  cells  used  in  the 
Denmark  test  facility  model.  Figure  C.15  is  an  example  of 
the  EARTH  output  for  the  Q1  code  in  Appendix  A. 

Program  run  time,  program  cost  and  EARTH  output  are 
directly  affected  by  the  number  of  cells  in  the  model  of 
the  physical  domain.  By  increasing  the  number  of  cells  in 
the  X,  y,  or  Z  direction  zones,  the  EARTH  program  requires 
more  "SLAB-WISE"  Computations  thereby  increasing  computer 
CPU  requirements  and  costs.  By  increasing  the  number  of 
cells  in  the  X  or  Y  directions  (i.e.,  NX  or  NY),  the  number 
of  cells  and  therefore  the  number  of  calculations  is 
increased.  For  example  for  a  fixed  NY,  SAY  16,  an  increase 
of  1  cell  in  the  X  direction  will  result  in  7  (number  of 
dependent  variables)  times  16  =  102  additional  equations 
(and  unknowns) .  An  increase  of  an  additional  cell  in  the  Z 


direction  will  result  in  7  times  NX  times  NY  additional 
equations.  For  example  for  NX=9  and  NY=16,  increasing  NZ 
from  28  to  29  will  result  in  an  additional  1,008  equations. 
Therefore,  deciding  to  use  a  “COURSE''  or  "FINE"  model  of 
the  physical  domain  is  determined  by  the  user's  time  and 
financial  constraints  and  the  desired  information  from  the 
flow  simulation  program.  The  Denmark  test  facility  model 
has  NX=9,  NY=16  and  NZ=28  for  a  grand  total  of  4,032  cells. 
Since  each  cell  has  7  dependent  variables,  28,224  linear 
algebraic  equations  must  be  solved  iteratively  to  achieve 
problem  convergence.  This  is  a  big  job  even  for  an  IBM 
mainframe  computer. 

For  the  purpose  of  this  research,  the  Denmark  test 
facility  was  modeled  with  a  "COARSE"  grid  with  the  intent 
of  achieving  problem  convergence  prior  to  using  a  "fine" 
grid.  Due  to  time  constaints  only  one  grid  was  used  in 
this  research.  Therefore,  problem  convergence  could  not  be 
identified  as  "GRID  INDEPENDENT". 

Although  the  Denmark  test  facility  model  grid  is 
considered  "COARSE",  many  of  the  cells  vary  in  size.  The 
number  and  size  of  cells  for  a  particular  region  of  the 
grid  was  chosen  based  upon  the  expected  activity  in  a 
region.  For  example,  at  the  primary  air  inlet,  one  would 
expect  primary  cooling  air  to  be  entering  at  a  relatively 
low  velocity  and  in  a  laminar  and  uniform  flow.  Then  as 
the  primary  air  approaches  the  jet  engine  intake,  fluid 
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velocities  should  be  increasing  in  magnitude  and  changing 
in  direction.  At  the  outlet  of  the  jet  engine  there  exists 
a  very  hot  and  a  very  fast  fluid  flow  mixing  with  entrained 
primary  cooling  air.  In  this  particular  region  with  a  lot 
or  activity,  a  less  ''COARSE”  grid  is  more  desirable  in 
order  to  incorporate  the  mixing  process  into  the  model .  As 
the  exhaust  gas  enters  the  augmentor  tube  it  mixes  with  the 
cooler  entrained  secondary  air.  Again,  a  less  "COARSE" 
grid  is  necessary  to  better  analyze  the  activity  involved. 

Figure  C.16  is  a  schematic  diagram  of  the  system  domain 
grid  for  the  X-Z  plane  while  Figure  C.17  represents  the  Y-Z 
plane.  The  two  figures  illustrate  a  less  "COARSE"  grid 
arrangement  at  the  jet  engine  exhaust,  augmentor  tube 
inlet,  and  secondary  air  inlet.  Again,  grid  refinement  for 
any  physical  domain  zone  is  accomplished  by  changing  one 
user  defined  variable  in  the  Q1  code. 

Although  the  Denmark  test  facility  was  modeled  as  a 
three-dimensional  structure.  The  body  fitted  coordinate 
option  of  PHOENICS  was  not  utilized.  As  such,  the 
cylindrical  jet  engine,  the  augmentor  tube,  and  the  exhaust 
stack  turning  vanes  were  modeled  as  rectangular  cubes.  The 
jet  engine  was  specified  as  a  solid  block  with  a  hot  mass 
source  at  one  end  and  a  mass  sink  at  the  other  end.  The 
jet  engine  was  modeled  with  CONPOR  and  PATCH-COVAL 
statements. 
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The  Denmark  test  facility  walls,  floor,  and  ceilings 
were  modeled  utilizing  CONPOR  and  PATCH-COVAL  statements. 
The  CONPOR  statements  were  given  a  zero  porosity  value 
thereby  preventing  mass  transfer  across  any  solid 
structure.  PATCH-COVAL  statements  were  necessary  for  the 
walls,  floor,  and  ceilings  in  order  to  introduce 
appropriate  friction  and  set  fluid  velocity  components  to 
zero  magnitude  at  the  structure  surfaces. 

B.  MODEL  BOUNDARY  CONDITIONS  AND  TURBULENCE  MODELS 

Boundary  conditions  for  the  Denmark  test  facility 
openings  to  the  atmosphere  were  specified  with  PATCH-COVAL 
statements.  Ambient  pressure  was  set  to  101,325  Pascals 
while  ambient  temperature  was  set  to  295  K. 

To  incorporate  turbulence  into  the  flow  simulation 
problem,  this  research  utilized  the  PHOENICS  K-E  turbulence 
model.  The  K-E  model  introduces  the  KE  (turbulence  kinetic 
energy)  and  EP  (kinetic  energy  dissipation  rate)  dependent 
variables  into  the  problem.  Turbulence  modeling  is  not  an 
exact  science  and  the  choice  of  the  K-E  model  for  this 
research  was  arbitrary.  The  PHOENICS  user's  manual 
(Ref.  5)  explains  in  detail  the  K-E  turbulence  model  and 
three  other  turbulence  models  available  to  the  PHOENICS 
user. 

Readers  may  refer  to  Appendix  A,  the  annotated  Q1  code, 
which  explains  in  detail  the  user  defined  variables  for  the 


Denmark  test  facility  model.  Columns  1  and  2  are  reserved 
for  the  beginning  of  executable  statements,  while  non¬ 
executable  statements  (comments)  begin  in  column  3  or 
beyond . 
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V.  ANALYSIS  OF  RESULTS 

A.  OBTAINING  A  FINAL  SOLUTION 

The  original  purpose  of  this  research  was  to  computer 
model  the  Royal  Danish  Air  Force  Gas  Turbine  Test  Facility 
using  the  PHOENICS  code.  The  computer  model  would  then  be 
used  to  compare  field  data  obtained  from  the  test  facility 
with  computer  generated  pressures,  temperatures,  and 
velocities.  If  the  computer  generated  output  compared  well 
with  the  field  data,  one  might  suggest  that  the  computer 
model  was  successful.  However,  if  field  data  from  a  test 
facility  was  non-existent  or  perhaps  very  limited  how  then 
could  the  success  of  a  computer  model  be  measured? 

The  field  data  from  the  Denmark  test  facility  was 
limited  in  scope.  Besides  the  mass  flow  described  in 
Chapter  II,  the  only  other  field  data  available  are  the 
test  bay  pressure  and  the  velocity  and  temperature  of  the 
engine  exhaust  at  one  point  in  the  augmentor  tube  and  at 
one  point  in  the  exhaust  stack.  Since  the  field  data  from 
the  Denmark  test  facility  was  limited,  a  comprehensive 
comparison  of  field  data  with  numerical  results  was  not 
possible.  To  investigate  the  EARTH  output  (numerical 
results)  and  conclude  that  the  numerical  results  were 
representative  of  the  test  facility,  this  study  examined 
the  EARTH  output  from  three  viewpoints:  common  sense 


30 


approach,  conservation  equation  balances,  and  problem 
convergence . 


B.  COMMON  SENSE  APPROACH 

The  Denmark  test  facility  was  modeled  with  five  major 
regions  or  zones.  Within  each  zone  a  common  sense  approach 
to  what  one  would  expect  the  gas  pressure,  temperature,  and 
velocity  conditions  to  be  was  taken.  The  five  major  model 
zones  are:  primary  air  inlet,  test  bay,  secondary  air 
inlet,  augmentor  tube,  and  exhaust  stack.  The  EARTH  output 
presented  in  Appendix  B  is  the  numerical  results  of  only  5 
slabs  in  the  major  model  zones. 

1.  Primary  Air  Inlet 

The  primary  air  inlet  is  comprised  of  the  following 
cells:  IX=1  to  IX=8,  IY=1  to  IY=12,  and  IZ=2.  Cool  ambient 
air  (295  K)  is  drawn  into  the  test  facility  through  the 
primary  air  inlet  (test  facility  doors) .  Pressure  drops  at 
the  test  facility  doors  were  very  small  with  a  maximum  drop 
of  about  12  Pascals.  Air  flow  through  the  primary  air 
inlet  was  predominantly  in  the  Z-direction  (W) .  Velocities 
ranged  from  1  to  2  m/sec  with  an  average  of  about  1.7 
m/sec.  The  horizontal  X-direction  velocity  (U)  and  the 
vertical  Y-direction  velocity  (V)  components  were 
negligible. 

With  an  average  gas  density  of  about  1.2  kg/m^  at 
the  test  facility  doors,  a  primary  air  mass  flow  rate  of 
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about  135  kg/sec  was  calculated.  Recall  from  chapter  IV 
that  the  Q1  program  modeled  only  one  half  of  the  test 
facility  due  to  the  Y-Z  plane  of  symmetry.  Therefore,  the 
135  kg/sec  primary  air  mass  flow  rate  is  for  one  half  of 
the  structure.  The  calculated  primary  air  mass  flow  rate 
for  the  entire  structure  would  be  about  270  kg/sec.  Of 
this  amount,  about  170  kg/sec  would  be  primary  cooling  air 
and  the  remaining  100  kg/sec  flows  through  the  engine. 

From  this  point  on,  references  to  mass  flow  rates  will  be 
for  the  entire  test  facility.  The  temperature  at  the  test 
facility  doors  remained  at  about  295  K.  At  the  test 
facility  doors,  one  would  not  expect  much  activity  and  a 
laminar  flow.  The  KE  and  EP  turbulence  dependant  variables 
were  negligible  which  indicated  laminar  flow. 

2 .  Test  Bay 

The  test  bay  comprises  cells  IX=1  to  IX=8,  IY=1  to 
IY=12,  and  IZ=2  to  IZ=9.  The  jet  engine  which  was  modeled 
as  a  2000  K,  100  kg/sec  sink-source  is  located  at  cells 
IX=1,  IY=5  to  IY=6,  and  IZ=6  to  IZ=7.  At  2000  K,  the 
density  of  the  engine  exhaust  gas  was  0.18  kg/m^ .  The  jet 
engine  provides  the  excitation  which  draws  ambient  air  in 
through  the  primary  air  inlet,  entrains  primary  cooling  air 
into  the  augmentor  tube  from  the  test  bay  and  entrains 


secondary  cooling  air  into  the  augmentor  tube  from  the 
secondary  air  inlets.  As  primary  air  moved  through  the 
test  bay  from  the  primary  air  inlet  to  the  sink  side  of  the 


jet  engine,  its  pressure  and  velocity  changed.  At  the  two 
cells  (IX=1,  IY=5  to  IY=6,  and  IZ=5)  before  the  engine 
inlet,  about  1.5  meters  upstream  of  the  inlet,  primary  air 
pressure  dropped  about  40  Pascals.  The  primary  air  W- 
velocity  component  was  about  3  m/sec  at  IZ=4,  about  3 
meters  upstream  of  the  engine  inlet.  The  horizontal  U- 
velocity  component  was  about  -4.5  m/sec  which  indicated 
that  air  was  being  drawn  into  the  engine  from  the  side 
walls  of  the  test  bay.  The  vertical  V-velocity  component 
was  about  -4  m/sec  above  the  engine  and  about  4  m/sec  below 
the  engine.  This  indicated  that  air  was  being  drawn  down 
from  the  test  bay  ceiling  into  the  engine  and  at  the  same 
time  being  drawn  up  from  the  test  bay  floor.  The  turbulent 
parameters  KE  and  EP  were  very  small  indicating  laminar 
behavior  at  the  engine  inlet.  The  temperature  at  the 
engine  inlet  remained  at  about  295  K. 

At  the  discharge  side  of  the  engine  (cells  IX=1, 

IY=5  to  IY=6,  and  IZ=8) ,  significant  activity  was  observed. 
The  pressure  increased  about  2100  Pascals.  In  adjacent 
cells  above/below  and  longside  the  engine  at  IZ==8,  the 
pressure  dropped  about  410  Pascals.  The  pressure 
difference  among  the  cells  caused  an  expanding  or 
diverging  effect  if  the  jet  engine  exhaust  gas.  The  U- 
velocity  component  of  the  gas  was  about  30  m/sec  while  the 
V-velocity  component  was  about  50  m/sec.  The  jet  engine 

exhaust  gas  expanded  out  towards  the  rear  wall  of  the  test 
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facility  as  opposed  to  going  directly  into  the  augmentor 
tube.  However,  due  to  the  large  magnitude  of  the  W- 
velocity  component,  about  430  m/sec,  the  resultant 
direction  of  the  engine  exhaust  was  primarily  into  the 
augmentor  tube.  Hot,  fast  moving  gas  at  the  engine  outlet 
mixed  with  the  cool,  slow  moving  entrained  primary  cooling 
air.  Due  to  this  mixing,  small  turbulence  resulted,  this 
premise  is  supported  by  the  magnitude  of  the  KE  and  EP 
parameters  in  the  region.  KE  had  a  maximum  value  of  about 
110  J/kg,  while  EP  had  a  maximum  value  of  about  5400 
J/kg/sec.  Due  to  this  mixing  of  engine  exhaust  and  primary 
cooling  air,  over  a  distance  of  about  1.5  meters,  the 
exhaust  gas  was  cooled  from  2000  K  to  about  1830  K  prior  to 
entering  the  augmentor  tube  (IZ=10) .  Over  the  same  1.5 
meter  distance,  the  exhaust  gas  W-velocity  component 
decreased  from  about  430  m/sec  to  about  400  m/sec. 

The  investigation  of  the  test  bay  was  localized  to 
a  small  region  around  the  jet  engine.  For  the  most  part 
the  activity  in  the  rest  of  the  test  bay  was  uneventful 
(i.e.  no  turbulence  and  insignificant  pressure, 
temperature,  or  velocity  changes) . 

The  above  observations  in  the  test  bay  and  at  the 
inlet  and  outlet  ends  of  the  jet  engine  are  realistic.  One 
would  expect  primary  air  velocity  components  and  pressure 
to  change  somewhat  prior  to  going  into  the  engine  and  more 
significantly  at  the  engine  outlet.  Additionally,  due  to 


the  primary  air  being  drawn  into  the  engine  a  low  pressure 
area  at  the  inlet  makes  good  sense.  While  at  the  outlet 
side  of  the  engine  a  high  pressure  is  expected.  At  the 
engine  outlet  one  would  expect  that  the  hot,  fast  engine 
exhaust  colliding  and  mixing  with  the  cool,  slow  entrained 
primary  cooling  air  would  result  in  some  turbulence. 

3 .  Secondary  Air  Inlet 

As  was  mentioned  earlier  in  chapter  II,  the 
secondary  cooling  air  plays  a  significant  role  in  the 
operation  of  a  jet  engine  test  facility.  The  secondary 
cooling  air  serves  the  purpose  of  providing  additional  cool 
(ambient)  air  for  mixing  with  the  hot  engine  exhaust.  For 
the  Denmark  test  facility  the  largest  mass  flow  rate  occurs 
through  the  secondary  air  inlets.  The  test  facility  data 
indicated  a  mass  flow  rate  of  about  340  kg/sec.  In  the 
structure  model,  secondary  cooling  air  enters  the  augmentor 
tube  at  cells  IX=5,  IY=1  to  IY=10,  and  IZ=10  to  IZ=12.  The 
secondary  cooling  air  had  an  average  U-velocity  component 
of  about  -14  m/sec,  which  indicated  that  secondary  cooling 
air  was  being  drawn  into  the  augmentor  tube  from  outside 
the  structure. 

A  maximum  pressure  drop  of  about  480  Pascals  at  the 
secondary  air  inlet  was  in  consonance  with  the  relatively 
high  U-velocity  component.  With  an  average  air  density  of 
about  1.0  kg/m^  at  the  secondary  air  inlet,  a  mass  flow 


rate  of  about  420  kg/sec  was  calculated.  Although  the 


secondary  air  inlet  had  a  significant  pressure  drop  and  a 
relatively  large  U-velocity  component  (about  -14  m/sec)  it 
did  not  demonstrate  turbulent  behavior  based  upon  the 
relatively  low  KE  and  EP  parameters.  The  temperature  at 
the  secondary  air  inlet  indicated  a  temperature  rise  of 
about  15  K.  The  results  showed  that  warm  air  from  the 
exhaust  stack  was  flowing  into  the  secondary  air  inlet 
region. 

With  a  primary  cooling  air  mass  flow  rate  of  about 
170  kg/sec  and  a  secondary  cooling  air  mass  flow  rate  of 
about  420  kg/sec,  an  Augmentation  Ratio  of  about  5.9:1  was 
calculated. 

4 .  Auamentor  Tube 

Three  regions  of  the  augmentor  tube  were  chosen  for 
investigation:  The  region  IZ=10  to  IZ=16,  the  slab  IZ=20, 
and  the  slab  IZ=22.  At  the  beginning  of  the  augmentor 
tube,  IZ=10  to  IZ=16,  a  large  mass  of  primary  and  secondary 
cooling  air  collided  and  mixed  with  the  hot,  fast  exhaust 
gas.  As  a  result,  significant  turbulence  occurred.  A 
maximum  KE  value  of  about  7,100  J/kg  was  observed  while 
maximum  EP  was  about  500,000  J/kg/sec.  In  the  8.5  meters 
that  the  engine  exhaust  traveled,  it  cooled  from  2000  K  to 
about  920  K.  The  engine  exhaust  W-velocity  component 
decreased  from  about  430  m/sec  to  about  125  m/sec.  The  U 
and  V-velocity  components  were  considered  negligible. 
Figures  C.18  and  C.19  are  PHOENICS  enthalpy  (Hi)  and 
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velocity  (Wl)  profiles  for  the  cells  IX=1,  IY=1  to  IY=10, 
and  IZ=16.  The  reader  is  reminded  at  this  point  that  an 
enthalpy  value  divided  by  a  specific  heat  value  of  about 
1000  J/(kg  K)  will  give  a  good  temperature  approximation 
within  a  cell.  The  profiles  indicated  minimum  temperature 
and  W-velocity  values  at  the  bottom  (IY=1)  and  top  (IY=10) 
of  the  augment or  tube  and  maximum  temperature  and  W- 
velocity  values  at  the  centerline  of  the  augmentor  tube. 
These  temperature  and  W-velocity  distributions  make  good 
sense.  Since  the  location  iy=5  is  at  the  height  of  the 
engine  outlet,  one  would  expect  the  engine  exhaust  to  be 
hotter  and  to  be  moving  faster  here  than  at  the  top/bottom 
of  the  augmentor  tube.  Figures  C.20  and  C.21  are 
temperature  and  W-velocity  profiles  for  the  cells  at  IX=1 
to  IX=5,  IY=5,  and  IZ=16.  Although  these  profiles  are  at 
the  same  slab  location  (IZ=16)  as  figures  C.18  and  C.19, 
the  examination  of  the  variables  is  from  the  center  of  the 
augmentor  tube  to  the  side  wall.  These  distributions 
indicated  maximum  values  of  the  variables  occurred  at  the 
center  of  the  augmentor  tube  and  minimum  values  at  the  side 
walls.  These  profiles  also  make  good  sense.  In  the  cells 
at  the  augmentor  wall  (IX=5,  IY=1  to  IY=10,  and  IZ=16) ,  the 
secondary  cooling  air  entering  the  augmentor  tube  has  only 
mixed  with  the  engine  exhaust  for  about  a  3  meter  distance 
down  the  augmentor  tube. 


At  augmentor  tube  location  IZ=20  (about  22  meters 
downstream  of  the  engine  outlet) ,  the  total  mass  flow  was 
examined.  The  W-velocity  component  ranged  from  about  ?5 
m/sec  to  75  m/sec  with  an  average  of  about  45  m/sec.  The  U 
and  V-velocity  components  were  considered  to  be  negligible. 
With  an  average  gas  density  of  about  0.6  kg/m^,  the  total 
mass  flow  rate  was  about  700  kg/sec.  At  augmentor  tube 
location  IZ=20,  there  was  stil?  some  turbulence,  which 
indicated  continued  mixing  of  the  primary/secondary  cooling 
air  and  the  engine  exhaust.  The  maximum  KE  value  was  abort 
740  J/kg  while  EP  maximum  was  about  10,400  J/kg/sec.  The 
temperature  ranged  from  about  520  K  to  about  635  K  with  an 
average  of  about  550  K.  Figures  C.22  and  C.23  are 
temperature  and  W-velocity  profiles  for  the  cells  at  IX=l, 
IY=1  to  IY=10,  and  IZ=20.  Figures  C.24  and  C.25  are 
temperature  and  W-velocity  profiles  for  the  cells  at  IX=1 
to  IX=5,  IY=5,  and  IZ=20.  A  more  uniform  distribution  of 
temperature  and  W-velocity  was  observed  across  the  slab  at 
IZ=20  than  at  slab  IZ=16.  This  is  reasonable  since  the 
total  mass  of  primary  and  secondary  cooling  air  and  engine 
exhaust  has  mixed  for  a  distance  of  about  22  meters 
downstream  of  the  engine  outlet. 

In  the  augmentor  tube/exhaust  stack  transition 
region  (IZ=22) ,  the  test  facility  field  data  indicated  the 
W-velocity  component  on  a  centerline  with  the  jet  engine 
was  about  60  m/sec  and  the  average  temperature  was  about 
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650  to  675  K.  The  numerical  results  indicated  a  W-velocity 
component  of  about  63  m/sec  and  a  temperature  of  about 
615  K.  In  a  distance  of  about  4  meters  from  IZ=20  to 
IZ=22,  turbulence  decreased  by  about  50  percent.  KE  was 
about  400  J/kg  while  EP  was  about  4,500  J/kg/sec. 

5 .  Exhaust  Stack 

Field  data  from  the  Denmark  test  facility  indicated 
a  temperature  rise  of  about  20  K  from  the  bottom  of  the 
exhaust  stack  to  the  area  under  the  exhaust  stack  baffle 
material.  In  the  numerical  model  a  temperature  drop  of 
about  50  K  was  observed.  The  Denmark  test  facility  was 
modeled  without  baffle  material  at  the  top  of  the  stack. 

Gas  velocities  at  the  top  of  the  computer  model  stack  were 
for  the  most  part  in  the  V-direction.  The  V-velocity 
ranged  from  about  10  to  45  m/sec  with  an  average  of  about 
20  m/sec.  Turbulence  at  the  top  of  the  exhaust  stack  was 
considered  to  be  negligible.  Average  temperature  at  the 
top  of  the  exhaust  stack  was  about  550  K.  Figures  C.26, 
C.27,  C.28,  and  C.29  are  temperature  and  V-velocity 
profiles  for  the  top  (IY=13)  of  the  exhaust  stack.  Figures 
C.26  and  C.27  are  for  the  cells  IX=1,  IY=13,  and  IZ=2l  to 
IZ=26  (along  the  Z  axis).  Figure  C.26  indicated  highest 
temperature  values  towards  the  rear  (IZ=25)  of  the  exhaust 
stack.  Additionally,  Figure  C.27  indicated  that  the 
V-velocity  fluctuated  (low  to  high  to  low  ...)  from  IZ=21 
to  IZ=26.  It  is  suggested  that  the  manner  in  which  the 


turning  vanes  were  computer  modeled  caused  a  lack  of 
uniformity  in  the  temperature  and  V-velocity  distribution 
for  this  stack  region. 

Figures  C.28  and  C.29  are  the  temperature  and  V- 
velocity  distributions  for  the  cells  at  IX=1  to  IX=5, 

IY=13,  and  IZ=23  (from  the  center  of  the  exhaust  stack  out 
to  the  side  wall) .  These  figures  indicated  fairly  uniform 
temperature  and  V-velocity  profiles  throughout  the  cells  in 
this  region. 

C.  COMPUTER  DATA  AND  FIELD  DATA  COMPARISON 

The  difference  between  field  data  results  and  numerical 
results  can  be  attributed  to  several  things;  baffle 
material  (not  included  in  the  numerical  model) ,  turning 
vanes  (crudely  modeled  as  rectangular  cubes) ,  turbulence 
model  (modeling  is  not  an  exact  science) ,  and  model  grid 
(grid  independence  not  verified) . 

The  computer  model  total  mass  flow  rate  of  about  700 
kg/sec  is  about  15  percent  larger  than  the  600  kg/sec  total 
mass  flow  rate  of  the  Denmark  test  facility.  It  is 
suggested  that  not  including  the  baffle  material,  in  the 
computer  model  primary  and  secondary  air  inlets  and  the 
exhaust  stack,  resulted  in  larger  numerical  mass  flow  rates 
than  the  test  facility  mass  flow  rates.  Test  facility 
baffle  information  was  not  available  for  this  research. 


Field  data  indicated  a  maximum  pressure  drop  in  the 
test  bay  of  about  500  Pascals.  The  computer  model 
calculated  a  maximum  pressure  drop  of  about  410  Pascals  in 
the  cells  adjacent  to  the  jet  engine  outlet. 

As  was  mentioned  earlier,  the  exhaust  gas  temperature 
and  W-velocity  component  in  the  augmentor  tube/exhaust 
stack  transition  region  compared  very  well  with  the  field 
data . 

D.  CONSERVATION  EQUATION  BALANCES 

The  PHOENICS  output  code  enables  the  investigator  to 
determine  the  degree  conservation  of  mass  has  been 
achieved.  The  PHOENICS  output  symbol  for  a  mass  flow 
source  is  "Rl".  By  an  algebraic  summing  of  the  mass  flow 
sources,  a  conservation  of  mass  flow  calculation  can  be 
made.  For  a  total  accounting  of  all  the  mass  flow  sources 
in  a  computer  model,  the  algebraic  sum  of  the  sources 
should  be  zero.  For  the  Denmark  test  facility  this  sum  was 
equal  to  0.02  kg/sec.  If  this  amount  is  compared  to  the 
100  kg/sec  mass  flow  through  the  jet  engine  there  is  a  0.02 
percent  error.  Compared  to  the  700  kg/sec  mass  flow 
calculated  for  the  augmentor  tube  the  percent  of  error  is 
reduced  to  0.003.  Based  on  this  analysis,  a  mass  balance 
or  a  continuity  of  flow  exists  in  the  Denmark  test  facility 


model.  A  mass  balance  or  a  continuity  of  flow  is  a 
necessazy  but  not  sufficient  condition  for  problem 
convergence . 
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E.  CONVERGENCE 

Chapter  III  described  the  iterative  process  which 
adjusts  the  values  of  the  dependent  variables  during  the 
solution  process.  During  the  process,  imbalances  occurred 
between  the  left  and  right  hand  side  of  the  dependent 
variable  discretized  algebraic  equations.  These  imbalances 
are  called  residuals.  Problem  convergence  would  be 
obtained  when  all  residuals  were  zero  (i.e.  the  dependent 
variable  discretized  algebraic  equations  were  in  balance) . 
The  residual  goes  to  zero  as  the  number  of  sweeps  goes  to 
infinity.  With  numerical  solutions,  users  must  employ  a 
variety  of  methods  to  determine  problem  convergence  . ince 
residuals  can  not  decrease  to  zero,  as  a  practical  matter. 

With  the  PHOENICS  output  (Appendix  B) ,  this  research 
utilized  a  local  and  a  global  analysis  of  the  dependent 
variables  to  determine  problem  convergence.  A  local  or 
"spot  value"  analysis  can  be  accomplished  with  the  PHOENICS 
output  by  designating  cells  to  be  monitored.  Since  the 
test  facility  field  data  referenced  the  exhaust  gas 
temperature  and  velocity  for  a  location  in  the  augmentor 
tube/exhaust  stack  transition  region,  the  computer  model 
spot  checked  the  dependent  variable  values  at  the  same 
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location  (the  cell  at  IX=1,  IY=5,  and  IX=22) .  From  the 
initial  program  run  to  the  final  sweep  (about  10,000),  the 
spot  values  were  closely  monitored  and  compared.  Local 
convergence  at  this  cell  occurred  when  little  or  no  change 
in  the  dependent  variable  values  was  observed  from  sweep  to 
sweep.  Local  or  spot  value  convergence  of  dependent 
variables  for  one  cell  is  a  necessary  but  not  sufficient 
condition  for  problem  convergence.  For  this  reason,  the 
results  at  "ALL"  cells  were  spot  checked  during  the 
solution  process. 

During  the  problem  solution  process .  the  PHOENICS  code 
sums  the  dependent  variable  whole-field  residuals.  An 
examination  of  the  whole-field  sums  provides  a  means  of 
global  investigation  to  deteirmine  problem  convergence.  If 
whole-field  residuals  decrease  towards  zero,  with 
increasing  iterations,  problem  convergence  is  suggested. 
However,  if  the  residuals  increase,  with  increasing 
iterations,  problem  divergence  is  occurring  and  some 
adjustment  of  the  PHOENICS  relaxation  parameters  is 
required.  Global  analysis  is  also  accomplished  by 
examination  of  the  energy  balance  of  the  system.  The 
dependent  energy  variable  enthalpy  (HI)  residual  is  first 
divided  by  1,000,000.  This  result  is  then  divided  by  the 
absolute  value  of  the  sum  of  the  enthalpy  (HI)  net  sources. 
If  the  quotient  decreases  to  less  than  about  0.05  (5 
percent  error) ,  problem  convergence  is  suggested.  If  the 


quotient  continues  to  increase,  problem  divergence  is 
occur ing  and  PHOENICS  relaxation  adjustment  is  again 
required. 

F.  SUMMARY 

Although  there  was  limited  test  facility  field  data  to 
compare  with  numerical  results,  it  is  suggested  that  the 
computer  modeling  of  the  Denmark  gas  turbine  test  facility 
with  PHOENICS  code  was  successful.  The  field  test  data 
compared  well  with  the  numerical  results.  During  the 
common  sense  examination,  the  conditions  that  existed  in 
the  five  major  regions  of  investigation  made  good  sense  and 
were  indicative  of  what  one  might  expect  in  a  gas  turbine 
test  facility.  The  total  mass  flow  continuity  check 
indicated  less  than  0.003  percent  error.  Finally,  local 
and  global  convergence  was  indicated.  A  system  energy 
balance  was  within  a  5  percent  error. 


VI.  RECOMMENDATIONS 


A.  INTRODUCTION 

It  is  believed  the  PHOENICS  code  is  a  useful  tool  for 
the  solution  of  gas  turbine  flow  simulation  problems. 
However,  since  the  Denmark  test  facility  has  not  "FINELY" 
modeled  the  following  computer  model  enhancements  are 
suggested:  BFC,  PHOTON,  alternate  turbulence  models,  grid 

refinement,  and  baffle  modeling. 

B.  BFC 

Body  fitted  coordinates  (BFC)  are  necessary  to 
incorporate  the  cylindrical  and  curved  surfaces  of  the  jet 
engine,  augmentor  tube,  and  turning  vanes.  With  the  proper 
geometric  shapes  in  the  test  facility  model,  mass  flow  in 
the  augmentor  tube  inlet,  the  secondary  air  inlets,  and  the 
exhaust  stack  will  be  more  representative  of  the  actual 
test  facility. 

C.  PHOTON 

The  implementation  of  the  PHOENICS  PHOTON  option  will 
greatly  enhance  numerical  results  investigation.  A  two  or 
three  dimensional  representation  of  a  dependent  variable 
field  is  much  more  useful  than  a  cell  by  cell  examination 
of  numbers. 
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D.  ALTERNATE  TURBULENCE  MODELS 

One  or  more  of  the  alternate  PHOENICS  turbulence  models 
should  be  implemented  and  investigated.  Since  it  was 
suggested  that  turbulence  modeling  is  not  an  exact  science, 
verifying  the  effect  of  the  turbulence  model  on  the 
solution  is  crucial .  It  may  be  that  the  turbulence  model 
selected  does  not  greatly  effect  the  solution. 

E .  GRID  REFINEMENT 

Increasing  the  number  of  cells  (grid  refinement)  in  the 
system  will  enhance  the  problem  solution  process.  During 
the  slab-sweep  iterative  process,  the  dependent  variable 
discretized  algebraic  equations  will  be  solved  in  a  "LESS 
ABRUPT"  manner  which  will  promote  problem  convergence. 

Grid  refinement  is  also  necessary  to  verify  the  solution  is 
grid  independent. 

F.  BAFFLE  MODELING 

Finally,  by  utilizing  additional  PHOENICS  CONPOR 
commands,  baffle  material  could  be  modeled  into  the 
primary/secondary  air  inlets  and  the  exhaust  stack.  By 
experimenting  with  different  values  of  porosity,  the  total 
mass  flow  rate  of  the  computer  model  could  more  closely 
represent  the  total  mass  flow  rate  of  the  actual  test 
facility. 
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APPENDIX  A 
Q1  CODE 


« 


TALK=F>RUN(1,1J 


DENMARK 

JET-ENGINE  TEST  FACILITY  COMPUTER  MODEL 

* 

NON-PHOENICS  USER-DEFINED  VARIABLES 

» 

XLBLl 

THE  FOLLOWING  VARIABLES  ARE  USED  TO  DESCRIBE  THE 

« 

XLMIDl 

PHYSICAL 

CHARACTERISTICS  OF  THE  DENMARK  TEST 

« 

XLM1D2 

FACILITY. 

THE  FIRST  LETTER  OF  THE  VARIABLE 

XLMID3 

SPECIFIES 

THE  DIRECTION  IN  THE  COORDINATE 

« 

XLBL2 

SYSTEM. 

XLBL3 

THE  FOLLOWING  CONVENTIONS  ARE  USED  TO  DESCRIBE 

XLEXTl 

A  PARTICULAR  ZONE  OF  THE  MODEL: 

« 

XLEXT2 

L 

LENGTH  IN  A  COORDINATE  DIRECTION 

XNCBLl 

MID 

MIDDLE  OF  SAY  THE  AUGMENTOR  TUBE 

» 

XNCMDl 

BL 

BOUNDARY  LAYER 

« 

XNCMD2 

EXT 

EXTERIOR  ZONE  OF  STRUCTURE 

XNCMD3 

NC 

NUMBER  OF  CELLS  IN  A  DIRECTION 

» 

X^JCBL2 

INT 

INTERNAL  ZONE  OF  STRUCTURE 

XNCBL3 

« 

XNCEXl 

EXAMPLE: 

XLBLl  THE  LENGTH  OF  THE  FIRST 

« 

XNCEX2 

ZONE  IN  THE  X  DIRECTION. 

« 

YLMID 

» 

YLBH 

YLBL2 

YNCEX2  THE  NUMBER  OF  CELLS  FOR  THE 

« 

YLEXTl 

SECOND  EXTERNAL  ZONE  IN  THE  Y 

« 

YLEXT2 

DIRECTION. 

# 

YLBL3 

« 

YLEXT3 

ZLINTS  THE  LENGTH  OF  THE  FIFTH 

« 

YLBL4 

INTERNAL  ZONE  IN  THE  Z 

« 

YNCMIO 

DIRECTION. 

» 

« 

YNCBLl 

« 

YNCBL2 

YNCEXl 

* 

YNCEX2 

« 

YNCBL3 

» 

YNCEX3 

« 

YNC3L4 

2LEXT1 

2LINT1 

ZLIMT2 

* 

» 

ZLINT3 

ZLINT4 

* 

ZLINT5 

» 

ZLINT6 

« 

ZUEXT2 

« 

ZNCEXl 

* 

ZNCINl 

* 

ZNCIN2 

# 

* 

ZMCIN3 

» 

# 

ZNCIN4 

* 

2NCIN5 

2NCIN6 

2NCEX2 


« 

« 

TMPJET  TEMPERATURE  OF  JET  EXHAUST  IK)  « 

CPJET  SPECIFIC  HEAT  OF  JET  EXHAUST  rj/KG»KJ  « 

PRSJET  PRESSURE  OF  JET  EXHAUST  IN/M**2)  « 

GASCON  GAS  CONSTANT  OF  JET  EXHAUST  1N*M/KG*K)  if 

MASJET  MASS  FLOW  RATE  OF  JET  EXHAUST  I  KG/SEC)  * 

VELJET  VELOCITY  OF  JET  EXHAUST  IM/SEC)  * 

RHOJET  DENSITY  OF  JET  EXHAUST  (KG/MiHf3)  * 

« 

TMPAMB  AMBIENT  AIR  TEMPERATURE  (K)  * 

PRSAMB  AMBIENT  AIR  PRESSURE  (N/M**2)  * 

CPAMS  SPECIFIC  HEAT  AMBIENT  AIR  I J/KG»K  )  * 

« 

« 

« 

« 


« 

if  PHOENICS  ALLOWS  45  USER  DEFINED  REAL-NUMBER  VARIABLES 
* 

REAL(XLBLl,XLMIDl,XLMl02,XLMID5>XL8L2,XLBLS>XLEXTi>XLEXT2  ) 

REAL! YLMI0,yLBLl,YLBL2,YLEXTl,yLEXT2,YLBL5>YLEXT5,YLBL4  » 

REAL! 2LEXT1 ,2LIMT1 ,2LINT2 ,ZLINT3 ,2LINT4 ,2LINT5 .ZLINTb ,2LEXT2  1 
REAL! IMP JET, CPJET, PRSJET, GASCON. MASJET, VELJET, RHOJET) 

REAL!  TMPAMB,  PRSAr^B,CPAN3,DELT,  FACT) 

REAL(KEINJ,EPINJ,KEINA,EPINA) 

if  PHOENICS  ALLOWS  45  USER  DEFINED  INTEGER  VARIABLES 

INTEGERIXNCBL1,XNCM01,XNCMD2,XNCM03,XNC8L2,XNCBL5,XNCEX1,XNCEX2  1 
INTEGER! YNCMIO , YNCBLl ,YNCBL2 .YNCEXl ,YNCEX2 ,YNCBL5 , YNCEX3 ,YNCBL4  ) 
INTEG£RlZNCEXl,ZNCINl,ZrJC:N2,ZNCIN5,2NCIN4,cNCIN5,2fJCIN6,Zr;CEX2  ) 


« 

* 


» 


#  * 


#  * 

if  INITIALIZATION  OF  USER-DEFINED  VARIABLES  * 


*  X  DIRECTION  LENGTHS  (METERS) 

» 

» 

XLBL1=0.S 

XLMID1=0.S 

XLMID2=0.5 

XLMID3=0.S 

XLBL2=0.S 

XLBL3=5.4 

XLEXT1=2.0 

« 

»  NUMBER  OF  CELLS  IN  X  DIRECTION 

« 

» 

XNC8L1=1 

XNCM01=1 

XNCMD2=1 

XNCMD3=1 

XNCBL2=1 

XNCBL2=3 

XNCEX1=1 

* 

« 

*  Y  DIRECTION  LENGTHS  (METERS) 

« 

YLeLl=2.0 

YLMID=1.0 

YLSL2=2.0 

YLEXT1=3.85 

YLEXT2S2.1S 

YLeL3=6.0 

* 

» 

«  NUMBER  OF  CELLS  IN  THE  Y  DIRECTION 

« 

* 

* 

YNCBL1=4 

YNCMia=2 

YNCBL2=4 

YNCEX1=2 

YNCEX2=1 

YNCBL3=3 

* 

« 

*  Z  DIRECTION  LENGTHS  (METERS) 

« 

» 

* 

ZLEXT1=2.0 

ZLINT1=23.S 

ZLINT2=3.2 

ZLINT3=1.5 

ZLINT4=:6.0 

ZLINr5=I6.3 

ZLINT<j=:8.4 

ZLEXT2=4.0 

* 

•  NUMBER  OF  CELLS  IN  THE  Z  DIRECTION 

* 

» 

ZNCEX1=1 

ZNCIN1=S 

ZNCIN2=1 

2NCIN3=2 

ZNCIN4=6 

ZHCINS=S 

2HCIN6=6 

ZNCEX2=2 


« 


TEXT( DENMARK  JET-ENGINE  TEST  FACILITY  MODEL) 


*  GROUP  2.  TRANSCIENCE»  TIME-STEP  SPECIFICATION  « 


«  * 

*  GROUP  3.  X-DIRECTION  GRID  SPECIFICATION  • 

*  » 


*  * 

«  NX  s  THE  NUMBER  OF  CELLS  IN  THE  X  DIRECTION  (INTEGER)  » 

«  XULAST  s  DISTANCE  IN  THE  X  DIRECTION  (METERS)  * 

»  * 

NX=XNCSL  1 +XNCMD 1  tXNCMD  2  tXNCMDS  ♦:<NCBL  2+XNCBL3  ♦XNCEXl 

XUUAST=XL3L1+XLMID1*XLMID2*XLMID$+XL&L2<-XLBL5+XLEXT1 

*  * 

•  THE  METHOD-OF-PAIRS  HAS  BEEN  USED  TO  DEFINE  THE  NUMBER  « 

»  AND  SIZE  OF  CELLS.  THE  XFRAC.YFRAC,  AND  ZFRAC  SPECIFY  THE  X.Y,» 

«  AN'D  Z  COORDINATES  RESPECTIVELY.  THE  ODD  INDEX  SPECIFIES  THE  * 

»  NUMBER  OF  CELLS  IN  A  DIRECTION  AND  THE  EVEN  INDEX  SPECIFIES  THE» 

•  SIZE  OF  EACH  CELL.  THE  NEGATIVE  QUANTITY  FOR  THE  FIRST  ODD  * 

•  INDEX  IN  EACH  GROUP  IS  THE  RESERVED  WORD  WHICH  SPECIFIES  » 

»  "METHOD-OF-PAIRS"  * 

»  * 


XFRAC( 1  )  =  -XNCBLli 
XFRAC( 3  )=XNCMDli 
XFRAC( 5  l=XNCMD2) 
XFRAC( 7)=XNCMD3( 
XFRAC( 9  l=XNCBL2; 
XFRAC( 11  )=XNCBL3; 


XFRAC( 2  )=XLBL1/(XNC8L1«XULAST I 
XFRAC(  4  )=XLMI01/)  X)-(CMDl»tXULAST  ) 
XFRAC( 6  )=XLMID2/( XNCMD2«XULAGT ) 
XFRAC( 8  )=XLMID3/( XNCMD3»XULAST  ) 
XFRAC( 10  )=XLBL2/(XNC3L2»XULAST  ) 
XFRAC(  12  )=XLBL3/(X))CBL3»XULAST  ) 
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XFRACI13)=XNCEX1} 


XFRACI 14  )=XLEXT1/IXNCEX1*XULAST ) 


GROUP  4.  Y-DIRECTION  GRID  SPECIFICATION 


*  NY  =  THE  NUMBER  OF  CELLS  IN  THE  Y  DIRECTION  (INTEGER) 

*  YVLAST  s  DISTANCE  IN  THE  Y  DIRECTION  (METERS) 

« 

NY=YNCMI0*YNCBL1+YNCBL2+YNCEX1*YNCEX2+YNCBL3 

YVLAST=YLMID+YLBL1*YLBL2*YLEXT1*YLEXT2*YLBL3 

« 

*  METHOD  OF  PAIRS 


YFRAC(1)=-YNCBL1( 
YFRAC(3)sYNCMI0» 
YFRAC( S  )=YNCBL2> 
YFRAC( 7 1=YNCEX1> 
YFRAC(  9  |sY);CEX2i 
YFRAC( 11  l=YNCBL3t 


YFRAC( 2  )=YLBLl/( YNCBL1*YVLAST ) 
YFRAC( 4  )sYLMID/( YNCMID»YVLAST ) 
YFRAC( 6  )=YLBL2/( YNCBL2»YVLAST ) 
YFRAC( 8  l=YLEXTl/( YNCEX1*YVLAST ) 
YFRAC( 10  )=YLEXT2/( YNCEX2»YVLAST ) 
YFRAC( 12  i=VL3L3/( YNCBL3*YVLAST ) 


)>»««■»»«»»»»»»»»»»<»->«»«»)»>»»»»» It  »»»»»»»*»)<  WOK  » 
OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOKKOKOOKIIKKKKKOOOOOOOOOOOOOO 
O  O 

*  GROUP  S.  2-DIRECTION  GRID  SPECIFICATION  * 

o  o 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO  KKKWKKKOOKKKKKKKKWK  OOOOOOOOOOOOOOOOOO 

o  o 

*  NZ  =  THE  NUMBER  OF  CELLS  IN  THE  2  DIRECTION  (INTEGER)  • 

»  2HLAST  =  DISTANCE  IN  THE  2  DIRECTION  (METERS)  * 


NZ=ZNCEXlt2NCINl+ZNCIN2+ZNCIN3tZNCIN4+ZNCIN5+ZNCIN6+ZNCEX2 

ZKLAST=2LEXTl*ZLINTl*ZLINT2+ZLINT3tZLlNT4*ZLINT5+ZLINT6t2LEXT2 


ZFRAC( 1 l=-ZNCEXli 
ZFRAC( 3  )=ZNCIN1) 
ZFRAC(5)  =  ZNCIN2, 
ZFRAC(  7)=:ZNCi;(3( 
ZFRAC( 9)=ZNCIN4} 
ZFRAC(  11  )=ZNCIN5-, 
ZFRAC(  13  )=2r.'CIN6j 
ZFRAC(lS)=ZNCEX2i 


ZFRAC( 2  )=ZLEXTl/( ZNCEX1*ZHLAST ) 
ZFRAC( 4  )=ZLINTl/( ZNCIN1*ZWLAST ) 
ZFRACl 6  )  =  ZLINT2/( ZNCINZoZWLAST ) 
ZFRACl 8  )=ZLINTI/( ZNCIN3*ZWLAST ) 
ZFRAC(  10  )=2LINT4/(Z1CIN4*ZHLAST  ) 
2FRAC( 12  )=ZLINT5/) ZNCINSkZWLAST ) 
ZFRAC( 14  I=ZLINT6/(ZNCIN6*ZWLAST ) 
2FRAC1 16  )=2LEXT2/( ZNCEXZkZWLAST ) 


OOOOOKKKKKKKKKKKKOOOOOOKOOKOOOOKOKKOKOKOOOOOOOOKKOKKKKKK-KOOOOOOOKKO* 
K*00»000000»**»OOOK»0«0»»»0«0»0-KOK  K  K  K  K  KKOOOO-OOOOOOKKOOOOOOOOOOOOOOOO 


GROUP  6.  BODY-FITTED  COORDINATES  OR  GRID  DIRECTION 


GRCXIP  7.  VARIABLES  STORED.  SOLVED  AND  NAMED 


» 

» 

* 

S0LVE(P1,H1,U1,V1,)>(1) 

S0LUTN(P1,Y,Y,Y,N,N,N) 

STORE) RHOI.TMPl) 

« 

»  * 

*  * 

*  * 

*  GROUP  8.  TERMS  (IN  DIFFERENTIAL  EQUATION)  AND  DEVICES  * 

»  * 

*  * 

*  DEFAULT  TAKEN  * 

«  * 

«  * 

»  * 

**».***********»**l*****»*****************'<HMHf************ 
» »!>»»»»»»»»»■» 

«  * 

*  GROUP  9.  PROPERTIES  OF  THE  MEDIUM  (OR  MEDIA)  * 

*  * 
»»»»♦*»»«»»»»»#»<»»»»»*♦»»*»»  »X»»XXX)»XX»XX»«»»»XWX  *♦»««  X  X  »  X  X 

»  « 


TMP1A=TINY 

TMPAMB=295. 

PR3AMB=0. 

CPA);B=1004. 

TMP1B=1.0/CPAMB 

TMP1=:GRN02 

TMPJET-CCOO. 

CPJET  =  i00-4. 

PRSJBT-0. 

GASCOr4=287. 

RH01B=1. /GASCON 
MASJET=100. 

PRES30=101S2S. 

RHOJET  =  ( PRESSOtPRSJET  )/( GASCON*TMPJET ) 

VELJET=MA3JET/RH0JET 

RH01=GRNCS 

ENUL:1.0E-C5 

DRH1DP=GRN05 

TURMCOIKEMODL  ) 


STORE (VIST) 

K£INJ=0 . 5»( 0 . 005*VELJET )**2 
EPINJ=0 . 09«KEINJ«»1 .5/10. 05*0 . 5 ) 

K£INA=0.S*( 0 . 005*1 . 0 )**2 
£PINA=0 . 09*K£INA**1 .5/10.1) 

*  « 


GROUP  10.  INTERPHASE-TRANSFER  PROCESSES  AND  PROPERTIES 


DEFAULT  TAKEN 


* 


* 

» 


GROUP  11.  INITIALIZATION  OF  VARIABLE  OR  POROSITY  FIELDS 


* 

« 

» 


*  « 

FIINITl HI )=CPAMB»TMPAMB 
FIINITl TMPl  )sTMPAMB 
FIINIT1DEN1)=1.0 
FIINIT(P1)=0.0 
FIINIT(U1  )=0.0 
FIINITl VI  )=0.0 
FIINIT(H1)=0.0 
FIINITl KE  l*K£INA 
FIINIT;  E? liEPINA 

*  * 

«  « 

*  » 

*  BLOCKING  OUT  THE  MASS  OF  THE  JET  ENGINE  * 

*  » 
CONPOR(0.0,CELL,1,1,S,6,6,7) 

*  * 

*  WALL, FLOOR,  AND  CEILING  BOUNDARY  CONDITIONS.  RESTRICTING  FLOW* 

*  ACROSS  A  BOUNDARY  BY  USING  CONPOR  VALUES  FROM  0.0  (NO  FLOW)  TO* 

»  1.0  (MAX  FLOW).  BOUNDARY  FRICTION  TO  FLOW  IS  ALSO  IMPLEMENTED* 

*  USING  A  NEGATIVE  BOUNDARY  COORDINATE.  THIS  CONVENTION  CAUSES  * 

*  FLOW  PARALLEL  TO  THAT  BOUNDARY  FACE  TO  HAVE  FRICTION.  BOUNDARY* 

»  CONDITIONS  WILL  BE  DESCRIBED  BY  COMMENT  AND  PLANE.  « 

*  * 

*  « 

*  TEST  FACILITY  FOUNDATION  AND  GROUND  OUTSIDE  STRUCTURE  * 

*  * 
CCNPORl 0 . 0 , SOUTH , 1 ,NX , 1 ,1 ,1 ,NZ ) 

PATCHl FLOR, SHALL, 1, NX, 1,1,1,N2, 1,1) 

COVALl FLOR,U1,GKND2,0.0 ) 

COVALl FLOR,W1,GRN02,0.0 ) 

COVALl FLOR,KE,GRND2,GRND2 ) 

COVALl FLOR,EP,GRND2,GRND2) 


JET-ENGINE  TEST  BAY 


» 


*  BAY  SIDE  WALL  (Y-Z) 

» 

CONPOR(0.0,WEST,9,9,1,12,2,9) 

« 

*  BAY  CEILING  (X-Z) 


CC3NPORtO.O,SOUTH,1.8,13,13,2,9) 

*  BAY  WALL  REAR  (X-Y)  2  SECTIONS 

* 

CONPORI 0.0,LOW,6,8,1,12,-10.10  J 
CONPORI 0.0,LOW,1,S,11,12,-10,10) 

» 

*  EXHAUST  DUCT  INLET  TUBE 
« 

CONPORI 0.0, NORTH, If  2, 3, 3, ID, 14) 
PATCHITUB1,NWALL,1,2,3,3,10,14,1,1) 
COVAL(TUB1,U1,GRND2,0.0) 
COVALITU81,W1,GRN02,0.0) 

COVALI TUBl.KE ,GRN02,GRND2 ) 
COVALITUB1,EP,GRN02,GRHD2  ) 

CCMPORI 0.0,SOUTH,1,2,8,8,10,14I 
PATCH(TU32,SWALL,1,2,8,3,10,14,1,1) 
COVALI TU32 , Ul, GRND2, 0. 0 ) 

COVALI TUB2 ,H1 ,GRMD2 ,0 . 0 ) 

COVAL I TUB2  ,KE ,GSN02 ,GRN02  ) 
CCVALlTUB2,EP,3RND2,3RND2  I 
CONPORI 0.0, WEST, 3, 3, 4, 7, 10, 14) 
PATCHITU83,WWALL,3,3,4,7,10,14,1,1) 
COVALI  TUBZ,'.'1,GR, '.’22,0. 0  ) 
COVALITU33,W1,GRM02,0.0  1 
COVALI  TUB3  ,KE  ,GRriD2  ,GRN02  ) 

COVALI  rU83,EP,GR)ID2,GRND2) 

CONPORI 0.0, LOW,  1,5, 1,3, 10, 10  » 

CONPORI 0.0,  LOW, 1,5, 8, 10, 10, 10) 
C0NPORI0.0,L0l'l,3,S,4.7,10,10) 

« 

*  ScCONDARY-AIR  INLET  WALL  IX-Y) 


CONPORI 0.0, LCW, 6, 7, 1,10, 12, 12  ) 

PATCHI SAIW,LWALL,6,7,1,10,12,12,1,1) 

COVALI SAIW,U1,GRHD2, 0.0  I 
COVALI SAIW , VI ,GRND2 ,0.0) 

COVALI SAIW,KE  ,GRND2,GRN02  ) 

COVA  L I  SAI W , E  P , GRNO  2 , GRND  2  ) 

* 

*  SECONDARY-AIR  INLET  AND  PARTIAL  EXHAUST  DUCT  CEILING  I X-Z) 

* 

CONPORI 0.0, SOUTH, 1,7, 11, 11, 10, 12) 
PATCHISAIWPE,SWALL,1,7,11,11,10,12,1,1) 

COVALI 3AIWPE  ,U1,GRND2,0.0  ) 

COVALI SAIKPE,W1,GRN02, 0.0  t 
COVALI SAIWPE,KE,GRND2,GRND2 ) 


« 

« 


« 


» 

« 


« 


« 


» 

» 


« 


« 

« 


» 


« 

» 


« 


COVA  L ( SAIMPE , E  P , GRNO  2 , 6RN0  2 1 

*  * 

*  EXHAUST  DUCT  CEILING  (X-Z)  » 

»  « 

CONPORI 0.0, SOUTH, 1,5, 11,11, 13, 20) 

PATCH! EDC, SHALL, 1,5, 11, 11, 13, 20, 1,1) 

COVALI EDC ,U1 ,GRND2 ,0 . 0 ) 

COVALI EDC ,H1 ,GRN02 ,0 . 0 ) 

COVALI EDC ,KE ,GRND2 ,GRND2 ) 

COVALI E0C,EP,6RN02,GRND2) 

»  * 

*  EXHAUST-DUCT  HALL  lY-Z)  * 

¥  » 

CONPORI0.0,HEST,6,6,l,10,13,20 ) 

PATCHIEDHXZ,HWALL,6,6,1,10,13,20,1,1) 

COVALI EDHXZ,V1 ,GRND2 ,0 . 0 ) 

COVALI EOHXZ, HI, GRN02, 0.0) 

COVALI EDHXZ,KE,GRND2,GRN02  ) 

COVAL I  E DHXZ , EP , GRNO 2 , GRND 2  ) 

*  K 

*  EXHAUST-STACK  HALLS  IX-Y)  TWO  SECTIONS  * 

»  * 

CONPORI 0 . 0 , HIGH ,1,5,11,13,20,20) 

PATCHI£SWXY1,HHALL,1,5,11,13,20,20,1,1) 

C0VALIESWXY1,U1,GRND2,0.0 I 
COVALI ESHXYl, VI, GRND2, 0. 0 ) 

C  0 VA  L ! £  SKXY 1,KE,GRND2.GRND2) 

COVA  L I ESHXY 1 , £  P , GRND  2 , GRNO  2 ) 

C0HPOR(0.0,LOII,l,5,l,13,27,27) 

PATCHIESWXY2,LWALL,1,S,1,13,27,27,1,1) 

COVALI ESWXY2,U1,GRN02,0.0 ) 

COVALI ESWXY2,V1,GRND2,0.0 ) 

COVALI ESHXY2,KE ,GRND2 ,GRN02  ) 

COVALI ESWXY2,EP,GRN02,GRND2 ) 

*  * 

*  EXHAUST-STACK  WALL  lY-Z)  * 

■*  * 

CONPORI 0.0, WEST, 6, 6, 1,13, 21, 26  I 
PATCHIESWXZ,WNALL,6,6,l,i3,21,26,l,l) 

COVALI  ESW;';Z,V1,GRN02,0.0  t 
COVALI ESHX2,Wi,GRND2,0.0 I 
COVALI  ES;'<XZ,KE  ,GRN02,GRND2  1 
C0VALIESWXZ,£P,GRND2,GRND2) 

»  * 

*  TURNING  VANES  * 

*  * 
CONPORI 0.0, SOUTH, 1,5, -7, 7, 21, 22) 

CONPORI  0 . 0  .SOUTH  ,l,S,-3,3,23,2<t) 

CONPORI 0.0, LOW, 1,5, 7, 10, -23, 23 ) 

CONPORI 0.0, LOW, 1,5, 3, 6, -25, 25) 

K  »*«*****»*****»«*********»* 

»  « 

»  GROUP  12.  CONVECTION  AND  DIFFUSION  ADJUSTMENTS  * 


DEFAULT  TAKEN 


« 

* 


* 

« 

* 


GROUP  IS.  BOUNDARY  CONDITIONS  AND  SPECIAL  SOURCES 


» 

*  JET-ENGINE  MASS  SOURCE 

* 

PATCH! JETIN, HIGH, 1.1,5, 6, S.S, 1,1) 

COVALI JETIN,P1,FIXFLU,-MASJET*1.  ) 

COVALC JETIN,Hl,ONLYMS,SAME  ) 

PATCH! JETOUT, LOW, 1,1, 5, 6, 8, 8, 1,1) 

COVAL! JETOUT,Pl,FIXFLU,MASJET»l. ) 

COVAL! JETOUT,Wl,ONLYMS,VELJET ) 

COVAL!  JETOUT, HI, ONLYMS,CPJET»TMP JET) 

COVAL! JETOUT ,KE .ONLYMS ,KEINJ  ) 

CO VA  L ! JE  TOUT , E  P , ONL  YMS , E  PINJ ) 

» 

»  ATMOSPHERIC  BOUNDARY  CONDITIONS 

*  LAST  CELL  SURFACE  IN  X  DIRECTION 

• 

PATCHIXSKY,EAST,NX,NX,1,NY,1,NZ,1,1) 

COVAL! XSXY, PI, 0. 1, 0. 0  ) 

COVALI XSKY, HI, ONLYMS, CPAMB»TMPAMB) 

COVAL! XSKY ,KE , ONLYMS ,KEINA ) 

COVALI XSKY ,EP .ONLYMS ,EPINA ) 

*  LAST  CELL  SURFACE  IN  THE  Y  DIRECTION 

* 

PATCH! YSKY .NORTH , 1 ,NX ,NY ,NY , 1 ,N2,1 ,1 ) 

COVAL!YSKY,P1,0.1,O.OI 

COVAL! YSKY , HI .ONLYMS ,CPAMB»TMPAMB ) 

COVALI YSKY, KE, ONLYMS, KEINA  ) 

COVAL ! YSKY , E  P , ONL YMS , EPINA ) 

» 

«  FIRST  CELL  SURFACE  IN  THE  Z  DIRECTION 
* 

PATCH!ZSKY1,L0W,1,NX,1,NY,1,1,1,1) 

COVALI2SKV1,P1,0.1,0.0) 

COVAL! 23KY1 ,H1 .ONLYMS , CP AMS*TMPAMB ) 

COVAL IZSKYl.KE, ONL YMS.KEINA ) 

COVALI ZSKYl.EP, ONLYMS, EPINA) 

« 

*  LAST  CELL  SURFACE  IN  THE  Z  DIRECTION 
» 

PATCH! ZSKY2, HIGH, 1, NX, 1, NY, NZ.NZ, 1,1) 

C0VAL!ZSKY2,P1,0.1,0.0I 

COVAL! ZSKYZ , HI .ONLYMS ,CPAM3*TMPAMB  ) 

COVALI 2SKY2,KE .ONLYMS , KEINA ) 


« 


» 

» 

* 

* 


m 

« 


* 

* 


* 


» 


* 


COVAL ( ZSKY2 ,EP ,ONLYMS ,EPINA  J 


*  * 

*  GROUP  14.  DOWNSTREAM  PRESSURE  FOR  PARAB=T  » 

«  * 

a  * 

*  DEFAULT  TAKEN  » 

«  « 

«  « 

it  GROUP  IS.  TERMINATION  OF  SWEEPS  * 

w  » 

RESTRTt ALL) 

FSWEEP=1 

LSWEEP=2 

»  » 

»  it 

*  * 

«  it 

*  GROUP  16.  TERMINATION  OF  ITERATIONS  * 

*  it 

»»»ititit»itit»»itltit«)<l<itiHtitititi(KWWItitMMItWitiHtiHtitiHtiHtit«WWKitiHHHtitXXX>tM)mMMititititit 
it  it 

*  DEFAULT  TAKEN  * 

it  if 

ititififitititiHtitifiHHt«ititititifiHtifititiHtiH(«itiHtitiHfitiHtit»iHtifitiHfiHtitiHHfifiHtitifiHfif#iHfitifiHf 
it»<<tiHtitititiHtititititiHtititititiHtiti»i»ititiHHt»itititi»iHHtititi»iH»iHtitiHtit»itititiHtit  X  X  M  K  M  M  Xititititit 

*  » 

if  GROUP  17.  UNDERRELAXATION  DEVICES  it 

*  * 

itit»»itit<titititititititiHtititititititititititiHtitiHt»i»it  X  X  X  Xitit»itititit»itit»»iHtititiHtititiHHtitX  X  X  Xitit 
it  it 

DELT=50  .  /I  NZitO .  ZitVELJET  ) 

FACT=0.25 

SHIFTED  TO  FACT-0,S  AT  1550  SWEEPS 
FACT=0.5 

RELAX1P1,LINRLX,0.3) 

RELAXlUl,FALSDT,FACTitOELT) 

RELAXl  Vl,FALSDT,FACTitOELT  1 
RELAX*  W1 ,  FALSOT  ,FACTitOELT  ) 

RELAXl  KE  ,  FALSDT,  FACTitOELT  ) 

RELAXl  EP  ,  FALSDT ,  FACTitDE LT  1 

REMOVED  THE  2.5  FROM  HI  RELAX  AT  5000  SWEEPS 
RELAXl  HI , FALSDT, FACT^DELTitl .  0  ) 


GROUP  18.  LIMITS  ON  VARIABLES  OR  INCREMENTS  TO  THEM 


»  * 
VARMIN(P1)=-2»PRESS0 
VARMIN( HI  )=CPAMB*TMPAM8 
VARMIN(U1)=-VELJET 
VARMINl VI  )  =  -VELJET 
VARMINIWl  )=-VnLJET 
VARMINl TMPl )=TMPAMB 
VARMINl RHOl )=0.1 
VARMAX(Pl)=2*PRESSO 
VARMAXI HI )=CPJET»TMPJET 
VARMA.X(U1)=VELJET 
VARMAXI VI )=VEL JET 
VARMAXI HI )=VELJET 
VARMAXI TMPl )=TMPJET 
VARMAXI RHOl )=2.0 

*  * 

«  » 

»  * 


*  « 

*  GROUP  19.  DATA  COMMUNICATED  BY  SATELLITE  TO  GROUND  » 

#  * 

».  ^  ^  W*.»WWV..<WWWWW.^W^WS4.*V.WV^ 

*  * 

*  DEFAULT  TAKEN  * 

*  » 

it  * 

GROUP  20.  PRELIMINARY  PRINTOUT  * 


it  it 

it  it 

it  DEFAULT  TAKEN  * 

it  * 

)<  KitilititiHtitititiHtitit  M  M  »  *  )titititititititiHHtititititiHt»»it 
ititT**it-jtit#itititit^itititiHtiHtitititititit*itititititiRitiHRititit-»itiHtitititiHtitiHt4Htit#ititit*iHtiHt*itit*it 


GROUP  21.  PRINTOUT  OF  VARIABLES  * 

if  * 

ifitititititititit»iHtitititititititititititititititititititititititiHHtitit)(  X  M  M  M  W  W  W  N  MititititiHHtititititiHtiHtititit 
it  it 

*  DEFAULT  TAKEN  » 

it  it 

itititiHtititititit»ititititititititititititiHtit  Xitit  ititititititititititititit M  X  M  itit  itiHHtititiHtitiHHtititititiHtititit 


it 

it 


GROUP  22.  SPOT-VALUE  PRINTOUT 


« 

it 


IXM0N=1 

IYMON=S 

IZMON=22 

TSTSWP=10 


PATCH(PL0TS,PR0FIL,1,1,S,5,1,NZ,1>1) 

PLOT( PLOTS, PI, -3000. 0,2500.0) 

PATCHl PL0T3,PR0FIL,1,X,S,S,1,NZ,1,1) 

PLOT) PLOTS , 01,-20.0, SO. 0 ) 

PATCH) PL0T4,PR0FIL,1,1,S,S,1,NZ, 1,1) 

PLOT) PL0T4 , VI ,-40.0,60.0) 

PATCH)PL0T2,PR0FIL,1,1,S,S,1,NZ,1,1) 

PLOT) PLOT2,W1,-1S.O,SOO.O ) 

PATCH) PL0T1,PR0FIL,1,1,S,S,1,NZ,1,1 ) 

PLOT) PLOTl,  HI, 275000. 0,2100000.0) 

PATCH) PL0T6,PR0FIL,1,1,1,10,20,20,1,1) 

PLOT) PL0T6, HI, 300000. 0,1000000.0) 

PATCH) PLOT7,PROFIL,1,1,1,10,20,20,1,1) 

PLOT)PLOT7,H1,20.0,17S.O  ) 

PATCH) PLOTS, PROFIL, 1,5,5,5,20,20,1,1) 

PLOT) PLOTS , HI  ,300000 . 0,1000000 . 0  ) 

PATCH) PLOT?, PROFIL, 1,5, 5, 5, 20, 20, 1,1) 

PLOT) PLOT?, Wl, 20.0, 175. 0  ) 

PATCH) PLOTIO,  PROFIL,  1,1, 1,10, 16, 16, 1, 1 ) 

PLOT) PLOTIO , HI ,300000 . 0,1000000.0) 

PATCH) PLOTll, PROFIL, 1,1, 1,10, 16, 16, 1,1) 

PLOT (PL0T11,W1, 20.0, 175. 0 ) 

PATCHl PLCT12, PROFIL  ,1,5, 5, 5, 16, 16, 1,1) 

PLOTl PL0T12,  HI,  300000. 0,1000000.0) 

PATCH) PL0T13, PROFIL, 1,5, 5, 5, 16, 16, 1,1) 

PLOTl  PLOT13,)Hl,20.0,17S.0  ) 

PATCHl PL0T14, PROFIL, 1,1, 13, 13, 21, 26, 1,1 ) 

PLOTl PLOT14,Hl,S00000.0,60000O.O ) 

PATCHl PL0T15,PR0FIL,1,1,13,13,21,26,1,1 ) 

PLOTl PLOT1S,V1,10.0,4S.O  ) 

PATCHl PLOT16, PROFIL, 1,5, 13, 13, 23, 23, 1,1) 

PLOTl PLOT16,H1,500000.0,600000.0 ) 

PATCHl PLOTl 7,  PROFIL,  1,5, 13, 13, 23, 23, 1,1 ) 

PLOT) PLOTl 7, VI, 10.0 ,45.0 ) 

NL'.'!CLS  =  10<‘NXPRI,'J 
NPLT=10 
ITABL=3 
IPROF=3 

LUPR3=6 

*  * 


»  * 

*,tlt*»*l>j»*IHHt**lt«IHHt»ltlt*»lt»*«tHHt«IK»»lt»l»IHt*««ltlt»IHtlt*l(  >  K*  Kll  K  »  It  lH>ltlt#JHK 


CO 


STOP 


APPENDIX  B 
EARTH  OUTPUT 
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APPENDIX  C 

FIGURES 


EXHAUST 


Figure  C.l.  Denmark  Jet-Engine  Test  Facility 


HIGH  < -  LOW 


GROUP 


TITLE 


RUN'  TITLE  AND  OTHER  PRELIMINARIES 
TRANSIENCE:  TIME-STEP  SPECIFICATIONS 
X-DIRECTION  GRID  SPECIFICATION 
Y-DIRECTION  GRID  SPECIFICATIONS 
2-DIRECTION  GRID  SPECIFICATIONS 
BODY-FITTED  COORDINATES  OR  GRID  DISTORTION 
VARIABLES  STORED,  SOLVED  AND  NAMED 
TERMS  (IN  DIFFERENTIAL  EQUATIONS)  AND  DEVICES 
PROPERTIES  OF  THE  MEDIUM  (OR  MEDIA) 
INTERPHASE-TRANSFER  PROCESSES  AND  PROPERTIES 
INITIALIZATION  OF  VARIABLE  OR  POROSITY  FIELDS 
CONVECTION  AND  DIFFUSION  ADJUSTMENTS 
BOUNDARY  CONDITIONS  AND  SPECIAL  SOURCES 
DOWNSTREAM  PRESSURE  FOR  PARAB=T 
TERMINATION  OF  SWEEPS 
TERMINATION  OF  ITERATIONS 
UNDERRELAXATION  DEVICES 

LIMITS  ON  VARIABLES  OR  INCREMENTS  TO  THEM 
DATA  COMMUNICATED  BY  SATELLITE  TO  GROUND 
PRELIMINARY  PRINTOUT 
PRINTOUT  OF  VARIABLES 
SPOT-VALUE  PRINTOUT 

23  FIELD  PRINTOUT  AND  PLOT  CONTROL 

24  DUMPS  FOR  RESTARTS 


Figure  C.4.  PHOEKTCS  Data  Hi-niip  ^R^ag^rie5 


FILE:  COPYQl  EXEC  A1 


&TRACE 

-START 

STYPE  WHAT  ARE  THE  FILENAME  AND  FILETYPE  OF  THE  Q1  FILE  YOU  WANT 
STYPE  TO  STORE  ON  MVS?  USE  ENTER  TO  EXIT. 
aREAD  VAR  SQl  a02 

SIF  /SQl  EQ  /  8EXIT 

STATE  aoi  a02  A 

aiF  8RC  EQ  0  SGOTO  -NEXT 

8TYPE  FILE  801  8Q2  A  DOES  NOT  EXIST.  USE  ENTER  TO  EXIT. 

8READ  VAR  8 
8IF  /a  EQ  /  8EXIT 
8G0T0  -START 
-NEXT 

STATE  COPYQl  JCL  A 

8IF  8RC  EQ  0  8G0T0  -COPY 

8TYPE  TEMPLATE  FILE  COPYQl  JCL  A  DOES  NOT  EXIST  ; 

8TYPE  SEE  LAB  MANAGER. 

8EXIT  -2 

-COPY 

AFPPnP  AFVTT  -1 

COFYFILE  COPYQl  JCL  A  TEMPI  JCL  A  (REPL 
8  =  aCONCAT  OF  /YOURQl FI L ENAME/  SQl  / 

8STACK  TOP 

8STACK  CLOCATE/YOURQIFILENAME/ 
aSTACK  CHANGE  8  1  1 
8STACK  CLOCATE/SYSUTl/ 

8STACK  GET  8QI  8Q2  A 
8STACK  TOP 
SSTACK  FILE 
XEDIT  TEMPI  JCL  A 
EXEC  SUBMIT  TEMPI  JCL  A 
CONWAIT 
DESBUF 
8EXIT 


FILE:  COPYQl  JCL  A1 


//CXHUSH  JOB  (A555,9999), 'COPYTMVS',CLASS=A 
//  EXEC  PGM=IEBGENER 
//SYSPRIfIT  DD  SYSOUT  =  A 
//SYSIN  DD  DUMMY 

//SYSUT2  DO  DSN=MSS.SA555.PH0ENIC3.SRC(Y0URQ1FILENAME) ,DISP=SHR 
//SYSUTl  DD  X 
/» 

// 


Figure  C.6  COPYQl  EXEC  And  JCL  Files 


FILE:  RUHSAT  EXEC 


A1 


&TRACE 

STYPE^WHAT  IS  THE  FILENAME  OF  THE  Q1  FILE  YOU  WANT  TO  RUN  (SATELITE) 
STYPE  USE  ENTER  TO  EXIT. 

&READ  VAR  8Q1 

aiF  /&QI  EO  /  8EXIT 
X  STATE  &Q1  DATA  A 
X  8IF  SRC  E(?  0  8G0T0  -NEXT 

X  STYPE  FILE  8Q1  DATA  A  DOES  NOT  EXIST.  USE  ENTER  TO  EXIT. 

X  8READ  VAR  8 

X  SIF  /8  EQ  /  8EXIT 

X  8GOTO  -START 
-NEXT 

STATE  RUN5AT  JCL  A 

8IF  8RC  EQ  0  8G0T0  -COPY 

8TYPE  TEMPLATE  FILE  RUNSAT  JCL  A  DOES  NOT  EXIST  ; 

8TYPE  SEE  LAB  MANAGER. 

8EXIT  -2 

-COPY 

SFPPnP  XFY7T  -1 

COPYFILE  RUNSAT  JCL  A  TEMP2  JCL  A  (REPL 
8  =  8C0NCAT  OF  /  YOURQl FI L ENAME  /  8Q1  / 

8STACK  TOP 
8STACK  CHANGE  8*1 
8STACK  FILE 
XEDIT  TEMP2  JCL  A 
EXEC  SUBMIT  TEMF2  JCL  A 
CONWAIT 
DESBUF 
8EXIT 


Figure  C.7  RUNSAT  EXEC  File 


FILE:  RUU5AT  JCL 


//RXHUSH  JOB  (4555,9999), »*CHAM',CLASS=A,REGION=102AK 
//G02  EXEC  PGM=SATLITP 

//5TEPLIB  DD  D5I1  =  MSS  .  PHOEMICS .  LOAD,  DISP  =  SHR 
//LIST!  DD  DSN=MSS. PHOENICS. LIST! .DIR, DISP=(OLD) 

//LI5T2  DD  DSM=MSS. PHOEMICS. LIST2. DIR. DISP=(OLD) 

//L1ST3  DD  DSM=MSS. PHOEMICS. LISTS. DIR, D1SP=(0LD) 

//LIST4  DD  DSN=MSS. PHOEMICS. LIST4. DIR, DISP=(OLD) 

//Ql  DD  DSM=MSS. S4555. PHOEMICS. SRCCYOURQl FI L EMAME) ,DISP=( OLD) 
//PHHELP  DD  DSM=MSS.PHOEtllCS.PHHELP.DATA,DISP=COLD) 

//SCRL  DD  DSM=MSS. PHOEMICS. SCRL,DISP=(MEH, PASS), 

//  DCB=(RECFM=F,LRECL=80,BLK3IZE=80) , 

//  SPACE=(TRK, ( 12,3) ) ,UMIT=SYSDA 

//COPY  DD  DSN=MSS. PHOEMICS. COPY, DISP=(NEW, PASS), 

//  DCB=(RECFM=F,LRECL=80, BLKSIZE=6320), 

//  SPACE=(TRK, ( 12,3)) ,UMIT=SYSDA 


//SCRL  DD  DSM=MSS. PHOEMICS. SCRL,DISP=(MEH, PASS), 

//  DCB=(RECFM=F,LRECL=80,BLK3IZE=80) , 

//  SPACE=(TRK, ( 12,3) ) ,UMIT=SYSDA 

//COPY  DD  DSN=MSS. PHOEMICS. COPY, DISP=(NEW, PASS), 

//  DCB=(RECFM=F,LRECL=80, BLKSIZE=6320), 

//  SPACE=(TRK, ( 12,3)) ,UMIT=SYSDA 

//FT09F001  DD  SYSOUT=x 
//FT06F0Q1  DD  SYSOUT=» 

//FT14F001  DD  SYSOUT=« 

//G512  DD  DSM^MSS. PHOENICS. G512,DISP=(0LD) 

//G513  DD  DSM=MSS. PHOEMICS. G513,DISP=(0LD) 

//G514  DU  DSM=MSS. PHOEMICS. G514,DISP=(0LD) 

//(3515  DD  D3M  =  MSS  .PHOEMICS. G515,  DISP  =  (OLD) 

//G516  DD  DSM=M3S. PHOEMICS. G5l6,DISP=(OLD) 

//G518  DD  DSN=MSS. PHOEMICS. G513,DISP=(0LD) 

//G526  DD  DSM=MS3. PHOEMICS. G526,DI5P=(0LD) 

//G520  DD  DSM=MSS. PHOENICS. G520,DISP=(OLD) 

//G521  DD  DSM=MSS. PHOEMICS. G521,DISP=(0LD) 

//G522  DD  DSM=MSS. PHOEMICS. G522,DISP=(0LD) 

//G523  DD  DSN=MSS. PHOEMICS. G525,DISP=(0LD) 

//G524  DD  D3N=MSS. PHOENICS. G524.DISP=(0LD) 

//G527  DD  DSM=MSS. PHOEMICS. G527,DI5P=(OLD) 

//G528  DD  D5M=MSS. PHOEMICS. G528,DISP=(0LD) 

//G529  DD  DSM=M5S. PHOEMICS. G529,DISP=(0LD) 

//G530  DD  DSr)  =  MSS. PHOEMICS. G530,DISP  =  (0LD) 

//G531  DD  D5N=M5S. PHOEMICS. G551.DISP=(0LD) 

//G735  DD  DSN=M5S. PHOEMICS. G735,DI3P=(0LD) 

//G736  DD  D3r)=:MS3. PHOEMICS. G756  ,DISP  =  (OLD) 

//G737  DD  DSM=M53. PHOEMICS. G737,DI5P=(OLD) 

//DFIO  DD  D5M=MSS.S4555. PHOEMICS. DFIO, 

//  DISP=SHR 

//DF12  DD  DSM=MSS.S4555. PHOEMICS. DF12,DISP=SHR 
// 


Figure  C.8.  RUNSAT  JCL  File 


a  c 


FILE!  RUNEAR  EXEC 


A1 


8TRACE 

EXEC  SUBMIT  RUNEAR  JCL  A 


FILE:  RUNEAR  JCL  A1 


//EXHUSH  JOB  (A555. 9999), 'CHAM', CLASSIC, REGI0N=2048K 
//*MAIN  ORG=NPGVMl .4555P,LINES=C12) 

//G02  EXEC  PGM=EARTHP 

//STEPLIB  DO  DSN=MSS.PHOENICS.LOAD,DISP=SHR 
//DSTL  DD  DUMMY 

//DF09  DD  DSN=MSS. 34555. PHOENICS.DF09,DISP=SHR 
//DFOl  DD  DSN  =  8«DF01,DISP=(NEll,PASS), 

//  DCB=(RECFM=VBS,LRECL=19000, BLKSIZE=19006 ), 

//  SPACE=(TRK,(25,3)),UNIT=SYSDA 

//DF02  DD  DSN  =  S8DF02,DI5P  =  (NEll,PASS), 

//  DCB=<RECFM=VBS,LRECL= 19000, BLKSI2E= 190 06 ) , 

//  SPACE=(TRK,(25,3) ),UNIT=SYSDA 

//DF03  DD  DSN  =  88DF03,DISP  =  (NEll,PASS) , 

//  DCB=(RECFM=VBS,LRECL=1 9000, BLKSI2E=1 9006), 

//  SPACE=(TRK,(25,3) ),UNIT=SYSDA 

//DF04  DD  DSN*88DF04,DISP=(NEM,PASS), 

//  DCB= (RECFM=VBS,LRECL= 19000, BLKSIZE* 19006), 

//  SPACE=(TRK, (25,3) ),UNIT=SYSDA 

//DF05  DD  D5N=88DF05,DISP=(NEW,PASS), 

//  DCB= (RECFM=VBS,LRECL= 19000, BLKSI2E= 19 006 ) , 

//  SPACE=(TRK,(25,3) ),UNIT=SYSDA 

//DF07  DD  D5N=a8DF07, DISP=(MEW,PASS) , 

//  DCB=(RECFM=VBS,LRECL=19000, BLKSIZE=19006) , 

//  SPACE=(TRK, (25,3) ) ,UNIT=5Y5DA 

//DF08  DD  DSN  =  aaDF08,DISP  =  (NEIJ,PASS), 

//  DCB=(RECFM=VBS, LRECL=1 900 0, BLKSIZE= 19006 ) , 

//  SPACE=(TRK, (25,3) ) ,UNIT=SYSDA 

//DF15  DD  DSN=aSDF15,DISP=(NEH,PA5S), 

//  DCB=(RECFM=VBS,LRECL=19000, BLKSIZE=19006 ) , 

//  SPACE=(TRK, (25,3) ) ,UNIT=SY5DA 

//DF16  DD  D3N=MSS.S4555.PH0ENICS.DF16,DISP=SHR 

//DFIO  DD  DSN=MSS. 34555. PHOENIC5.DF10,DISP=(OLD) 
//DF12  DD  D3N=M5S, 34555. PH0ENIC3.DF12,DI3P=(0LD) 
//FT14F001  DD  3Y50UT=K 
//FT06F001  DD  SYSOUT=» 

// 


Figure  C.9  RUNEAR  EXEC  And  JCL  Files 


Figure  C.ll.  Test  Facility  Top  View  (X-Z  Plane) 


METERS 


X  GROUP  3.  X-DIRECTION  GRID  SPECIFICATION 
X  NX  =  THE  NUMBER  OF  CELLS  IN  THE  X  DIRECTION  (INTEGER) 
X  XULAST  *  DISTANCE  IN  THE  X  DOMAIN  (METERS) 


X 

NX=XNCBL1+XNCMD1+XNCMD2+XNCMD3+XNCBL2+XNCBL3+XNCEX1 

XULA3T=XLBL1+XLMID1+XLMID2+XLM1D3+XLBL2+XLBL3+XLEXT1 


XFRAC(1)=-XNCBL1; 

XFRAC(5)=XNCMD1; 

XFRAC(5)=XNCMD2; 

XFRAC:? ;=XNCMD3; 

XFRAC(9)=XNCBL2; 

XFRAC(11)=XNCBL3; 

XFRAC(13)=XNCEX1; 


XFRAC(2)=XLBL1/(XNCBL1XXULAST) 

XFRAC(4)=XLMID1/(XNCMD1XXULAST) 

XFRAC(6)=XLMID2/(XNCMD2XXULAST) 

XFRAC(8)=XLMID3/(XNCMD3XXULAST) 

XFRAC(10)=XLBL2/(XNCBL2XXULAST) 

XFRAC(12)=XLBL3/(XNCBL3XXULAST) 

XFRAC{14)=XLEXT1/(XNCEX1XXULAST) 


X 

X  GROUP  4.  Y-DIRECTION  GRID  SPECIFICATION 
X  NY  =  THE  NUMBER  OF  CELLS  IN  THE  Y  DIRECTION  (INTEGER) 
X  YVLAST  =  DISTANCE  IN  THE  Y  DOMAIN  (METERS) 

X 

NY=YNCMID+YNCBL1+YNCBL2+YNCEX1+YNCEX2+YNCBLS 

YVLAST=YLMID+YLBL1+YLBL2+YLEXT1+YLEXT2+YLBL3 


if 


YFRAC(l)  =  -YtlCBLl; 
YFRAC(3)=YNCMID; 
YFRAC(5)=YNCBL2; 
YFRAC(7)=YNCEX1; 
YFRAC(9)=YNCEX2; 
YFRAC(11)=YNCBL3; 
X 


YFRAC(2)=YLBL1/(YNCBL1XYVLAST) 

YFRAC(4)=YLMID/(YNCMIDXYVLAST) 

YFRAC(6)=YLBL2/(YNCBL2XYVLAST) 

YFRAC(8)=YLEXT1/(YNCEX1XYVLAST) 

YFRAC(10)=YLEXT2/(YNCEX2XYVLAST) 

YFRAC(12)=YLBL3/(YNCBL3xYVLA5T) 


X  GROUP  5.  2-DIRECTION  GRID  SPECIFICATION 
X  NZ  =  THE  NUMBER  OF  CELLS  IN  THE  Z  DIRECTION  (INTEGER) 
X  ZWLAST  =  DISTANCE  IN  THE  Z  DOMAIN  (METERS) 


X 

NZ=ZMCEX1+ZNCIN1+ZNCIN2+2NCIN3+ZNCIN4+ZNCIN5+ZNCIN6+ZNCEX2 

2WuA5T=ZLEXTl+ZLINTl+ZLI»T2+ZLINT3+ZLlNT4+ZLINT5+2LINT6+ZLEXT2 


X 

ZFRACd  )=-ZNCEXl; 

ZFRAC(3)=2NCIN1; 

ZFRAC(5)=2NCIN2, 

ZFRAC(7)=2NCIN3; 

ZFRACC9)=ZNCIN4; 

ZFRAC(11)=ZMCIN5; 

ZFRACC 13)=ZNCIN6; 

ZFRAC(15)=ZNCEX2; 


ZFRAC(2)=ZLEXT1/(2NCEX1XZWLA5T) 

ZFRAC{4)=ZLINTl/(ZNCINlxZWLAST) 

ZFRAC(6)=ZLINT2/(2NCIM2XZWLAST) 

ZFRAC(8)=ZLINT3/(ZNCIN3XZWLAST) 

ZFRAC(  10)=ZLINT4/(ZrJCIN4xZWLA5T) 

ZFRAC(12)=ZLINT5/(ZNCIN5XZWLAST) 

ZFRAC(14)=ZLINT6/(ZNCIN6XZWLAST) 

ZFRAC(16)=ZLEXT2/(ZNC£X2xZWLA5T) 


Figure  C.IA.  Q1  MEthod-Of-Pairs  Coding  Sequence 


GROUP  3.  X-DIRECTION  GRID  SPACING 
CARTES  =  T 
NX  =  9 

XULAST  =  9.900E+00 

METHOD  OF  PAIRS  USED  FOR  GRID  SETTING 
XFRAC  C  1)  =  -l.OOOE+00  ;XFRAC  (  2)  =  5.0S1E-02 

XFRAC  (  3)  =  l.OOOe+00  ;XFRAC  <  A)  =  S.051E-02 

XFRAC  (  5)  =  l.OOOE+00  ;XFRAC  (  6)  =  5.051E-02 

XFRAC  (  7)  =  l.OOOE+QO  ;XFRAC  <  8)  =  5.051E-02 

XFRAC  (  9)  =  l.OQOE+00  ;XFRAC  (  10)  =  S.051E-02 

XFRAC  (  11)  =  3.000E+00  ;XFRAC  C  12)  =  1.818E-01 

XFRAC  (  13)  =  l.OOOE+00  ;XFRAC  (  14)  =  2.020E-01 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
GROUP  4.  Y-DIRECTION  GRID  SPACING 
NY  =  16 

YVLAST  =  1.700E+01 


METHOD  OF  PAIRS  USED  FOR  GRID  SETTING 


YFRAC 

( 

1) 

S 

-4. OOOE+00 

; YFRAC 

( 

2) 

S 

2.941E-02 

YFRAC 

( 

3) 

S 

2.000E+00 

; YFRAC 

( 

4) 

s 

2.941E-02 

YFRAC 

( 

5) 

s 

4.000E+00 

;YFRAC 

< 

6) 

s 

2.941E-02 

YFRAC 

( 

7) 

s 

2.  OOOE+QO 

; YFRAC 

( 

8) 

s 

1.132E-01 

i'FRAC 

( 

9) 

s 

l.OQOE+QO 

; YFRAC 

( 

10) 

s 

1 .265E-01 

YFRAC 

( 

11) 

5 

3.0006+00 

; YFRAC 

( 

12) 

S 

1.176E-01 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
GROUP  5.  Z-DIRECTION  GRID  SPACING 

parab  =  F 
NZ  *  28 

ZWLAST  =  6.490E+01 

METHOD  OF  PAIRS  USED  FOR  GRID  SETTING 
ZFRAC  <  1)  =  -l.OOOE+00  ;ZFRAC  (  2)  =  3.082E-02 

ZFRAC  (  3)  =  5.000E+00  ;ZFRAC  (  4)  =  7.242E-02 

ZFRAC  (  5)  =  l.OOOE+OO  ;2FRAC  (  6)  =  4.931E-02 

ZFRAC  (  7)  =  2.000E+00  ;ZFRAC  (  8)  =  1.156E-02 

ZFRAC  (  9)  =  6.000E+00  ;2FRAC  (  10)  =  1.541E-02 

ZFRAC  (  11)  =  5.000E+00  ;ZFRAC  (  12)  =  5.023E-02 

ZFRAC  (  13)  =  6.000E+00  ;ZFRAC  (  14)  =  2.157E-02 

ZFRAC  (  15)  =  2.000E+00  ;ZFRAC  (  16)  =  3.032E-02 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


PATCHCPLOTIO  ,PROFIL,  1,  1,  1,  10,  16,  16,  1,  1) 

PLOTCPLOTIO  ,H1  ,  3.000E+05.  l.OOOE+06) 

TABULATION  OF  ABSCISSA  AND  ORDINATES... 

Y  HI 

2.500E-01  3.45AE+05 

7.500E-01  4.154E+05 

1.250E+00  5.411E+05 

1.750E+00  7.754E+05 

2.250E+Q0  9.650E+05 

2.750E+00  9.S99E+05 

3.250E+00  7.667E+05 

3.750E+00  5.340E+05 

4.250E+00  4.097E+05 

4.750E+00  3.613E+0S 


VARIABLE  HI 

MINVAL=  3.000E+05 
MAXVAL=  l.OOOE+06 

CELLAV=  6.074E+05 

1.00  .  + - + 

0.90  +  + 

0.80  +  X  \  + 

0.70  +  A  + 

0.60  +  X  + 

0.50  +  X  + 

0.40+  + 

0.30+  H  H  + 

0.20+  >s.  + 

0.10  H^X^ 

0  .1  .2  .3  .4  .5  .6  .7  .8  .9  1.0 

THE  ABSCISSA  IS  Y  MIM=  2.50E-01  MAX=  4.75E+00 


Figure  C.18.  H  vs  Y-Axis  at  IY=1  to  IY=10 
For  IX=1  And  IZ=16 


PATCHCPLOTll  ,PROFIL,  1,  1,  1,  10, 

PLOT(PLOTll  ,W1  ,  2.000E+01,  1.750E+02) 

TABULATION  OF  ABSCISSA  AND  ORDINATES... 


2.500E-01 
7.500E-01 
1 .250E+00 
1 .750E+00 
2.250E+00 
2.750E+00 
3.250E+00 
3.750E+00 
4.250E+00 
4 .750E+00 
VARIABLE 
MINVAL= 
MAXVAL= 
CELLAV= 


W1 

2.923E+01 
4.467E+01 
7 . 167E+01 
1.226E+02 
1 .696E+02 
1 .684E+02 
1 .203E+02 
6 .873E+01 
4.100E+01 
2.919E+01 
W1 

2.000E+01 
1 .75QE+02 
8.655E+01 


16, 


16, 


1) 


Figure  C.19.  W  vs  Y-Axis  At  IY=1  To  IY=10 
For  IX=1  And  IZ=16 


1  LG 


PATCH(PL0T12 

,PR0FIL, 

1,  5,  5, 

PL0T(PL0T12 

,H1  ,  3. 

OOOE+05,  l.OOOE+06 

TABULATION  OF  ABSCISSA 

AND  ORDINATES. . . 

X 

HI 

2.500E-01 

9.650E+05 

7 .500E-01 

7 .746E+05 

1 .250E+00 

5.78AE+05 

1 .750E+00 

A.848E+05 

2.250E+00 

4.530E+05 

VARIABLE 

HI 

MINVAL= 

3. OOOE+05 

MAXVAL= 

1 . 000E+Q6 

CELLAV= 

6 .512E+05 

1) 


0.00  _ 

0  .1  .2  .3  .5  .6  .7  .8  .9  1. 

THE  ABSCISSA  IS  X  .  MIN=  2.50E-01  MAX=  2.25E+00 


Figure  C.20.  H  vs  X-Axis  At  IX=1  To  IX=5 
For  IY=5  And  IZ=16 


1  '  1 


+  -f  O 


PATCH(PL0T13  ,PROFIL,  1,  5,  5.  5,  16,  16, 

PL0T(PL0T13  ,W1  ,  2.000E+01,  1.750E+02) 


TABULATION 

OF  ABSCISSA  AND  ORDINATES 

X 

W1 

2.500E-01 

1 .696E+02 

7.500E-01 

1 . 196E+02 

1 .250E+00 

6 .867E+01 

1 .750E+00 

3.822E+01 

2.250E+00 

2.315E+01 

VARIABLE 

W1 

MINVAL= 

2. OOOE+01 

MAXVAL= 

1 .750E+02 

CELLAV= 

8.384E+01 

1) 


Figure  C.21,  W  vs  X-Axis  At  IX=1  To  IX=5 
For  IY=5  And  IZ=16 


1  :;g 


PATCH(PL0T6  ,PROPIL,  1,  1,  1,  10,  20,  20,  1,  1) 

PL0T(PL0T6  .HI  ,  3.000E+05,  l.OOOE+06) 

TABULATION  OF  ABSCISSA  AND  ORDINATES... 

Y  HI 

2.500E-01  5.117E+05 

7.500E-01  5.S3AE+05 

1.250E+00  5.850E+05 

1.750E+00  6.171E+05 

2.250E+00  6.376E+05 

2.750E+00  6.331E+05 

3.250E+00  6.057E+05 

3.750E+00  5.715E+05 

4.250E  +  00  5.<i47E+05 

4.750E+00  5.315E+0S 


VARIABLE  HI 

MINVAL=  3.000E+05 
MAXVAL=  l.OOOE+06 
CELLAV=  5.791E+05 

1.00 

0.90  +  + 

0.80+  + 

0.70  +  + 

0.60  +  + 

0.50  +  _ H - H - H - -  + 

0.40  +  - - -  H  + 

0.30  KT — — H 

0 . 20  -T  + 

0.10+  + 

0.00 

0  .1  .2  .3  .4  .5  .6  .7  .8  .9  1.0 

THE  ABSCISSA  IS  Y  .  MIN=  2.50E-01  MAX=  4.75E+00 


Figure  C.22.  H  vs  Y-Axis  At  IY=1  To  IY=10 
For  IX=1  And  IZ=20 


PATCHCPL0T7  .PROFIL,  1,  1,  1,  10, 

PL0T(PL0T7  ,W1  ,  2.000E+01,  1.750E+02) 

TABULATION  OF  ABSCISSA  AND  ORDINATES... 


20,  20, 


1, 


I) 


2.500E-01 
7.500E-01 
1 .250E+00 
1 .750E+00 
2.250E+00 
2.750E+00 
3.250E+00 
3.750E+Q0 
4.250E+00 
4.750E+00 
VARIABLE 
MINVAL= 
MAXVAL= 
CELLAV= 


HI 

4.170E+01 
5.000E+01 
5.790E+01 
6 .663E+01 
7 .281E+01 
7 .200E+01 
6 .108E+01 
5.2g7E+01 
4.782E+01 
4 .620E  +  01 
HI 

2.  OOOE+01 
1 .750E+02 
5.691E+01 


1.00 

0.90 

0.80 

0.70 


0.60 
0 . 50 


0  .1  .2  .3  .4  .5  .6  .7  .8  .9  1.0 

THE  ABSCISSA  IS  Y  .  MIN=  2.50E-01  MAX=  4,75E+00 


Figure  C.23.  W  vs  Y-Axis  At  IY=l  To  IY=10 
For  IX=1  And  IZ=20 

1  ‘t 


PATCH(PL0T8  ,PROFIL,  1,  5.  5.  5. 

PLOTCPLQTa  ,H1  ,  3.000E+05,  l.OOOE+06) 

TABULATION  OF  ABSCISSA  AND  ORDINATES... 


20,  20, 


1, 


1) 


X  HI 

2.500E-01  6.376E+05 

7.500E-01  6.165E+05 

1.250E+Q0  5.890E+05 

1.750E+00  5.679E+05 

2.250E+00  5.577E+05 

VARIABLE  HI 

MINVAL=  3.000E+05 
MAXVAL=  l.OOOE+06 
CELLAV=  5.937E+05 


00 

90 

80 


0.70 
0.60 
0.50 
0.40 
0.30 
0.20 
0.10 
0 . 00 

^  0  -2  .3  .4  .5  .6  .7  89i 
THE  ABSCISSA  IS  X  .  MIN=  2.50E-01  MAX=  2.25E+0. 


.  +  ....  + _ + _ + _ + _ +  , 


.  + 

+ 

+ 

+ 

+ 

+ 
+ 
,  + 

.  0 


Figure  C.24.  H  vs  X-Axis  At  1X=I  TO  1X=5 
For  IY=5  And  IZ=20 

1 1 1 


PATCH(PL0T9  .PROFIL,  1,  5,  5,  5,  20,  20, 

PL0T(PL0T9  ,W1  ,  2.000E+01,  1.750E+02) 


TABULATION 

OF  ABSCISSA  AND  ORDINATES 

X 

W1 

2.500E-01 

7 .281E+01 

7 .500E-01 

6 .47AE+01 

1 .250E+00 

5.391E+01 

1 .750E+00 

4.524E+01 

2.250E+00 

4 . 084E+01 

VARIABLE 

W1 

MirJVAL  = 

2.000E+01 

MAXVAL= 

1 .750E+02 

CELLAV= 

5 . 551E+01 

1, 


1  .00 

+  .  . 

. + . . . 

.  .  .  + 

0 .90 

+ 

+ 

0.80 

+ 

+ 

0.70 

+ 

+ 

0.60 

0 . 50 

+ 

+ 

+ 

+ 

0 . 40 
0 . 30 

+ 

V^ — 

— 

+ 

+ 

0.20 

0,10 

+ 

+ 

M 

u_ 

0  ,00 

+  .  , 

+  . 

.  .  .  + 

0 

.  1 

.2 

.3 

.4 

.5 

.6 

.7 

.8 

.9 

1  .0 

THE  ABSCISSA  IS  X  MIM=  2.50E-01  MAX=  2.25E+00 


Figure  C.25.  W  vs  X-Axis  At  IX=1  To  IX=5 
For  IY=5  And  IZ=20 


n 


1  1  2 


PATCH(PL0T14  ,PROFIL,  1,  1.  13.  13, 

PL0T(PL0T14  ,H1  ,  5.000E+05,  6.000E+05) 

TABULATION  OF  ABSCISSA  AND  ORDINATES... 


21,  26, 


1, 


1) 


1  .  OOOE+00 
4.350E+00 
9.050E+00 
1 .375E+01 
1  .845E+01 
2 . 315E+01 
VARIABLE 
MINVAL= 
MAXVAL= 
CELLAV= 


HI 

5.429E+0S 

5.535E+05 

5.653E+05 

5.895E+05 

5.905E+05 

5.631E+05 

HI 

5 .  OOOE+05 
6  .  OOOE+05 
5 .633E+05 


0  .1  .2  .3  .4  .5  .6  .7  .8  .9  1.0 

THE  ABSCISSA  IS  2  .  MIN=  l.OOE+00  MAX=  2.31E+01 


Figure  C.26.  H  vs  Z-Axis  At  IZ=21  To  1Z=26 
For  IX=1  And  IY=13 


PATCH(PL0T15  ,PROFIL,  1,  1,  13,  13,  21,  26,  1, 

PLOTCPLOTIS  ,V1  ,  l.OOOE+01,  A.SOOE+Ol) 

TABULATION  OF  ABSCISSA  AND  ORDINATES... 

Z  VI 

l.OOOE+00  1.286E+01 


4.350E+00 

2.973E+01 

9 . 050E+00 

2.223E+01 

1 .375E+01 

3.062E+01 

I  .845E+01 

2.919E+01 

2.315E+01 

2.821E+01 

VARIABLE 

VI 

MINVAL= 

1 .  OOOE+01 

MAXVAL= 

4.500E+01 

CELLAV= 

2.547E+01 

1.00  + - +  . 

...+....+... 

.  +  . 

.  .  .  +  . 

+ 

0.90  + 

+ 

0.80  + 

+ 

0.7C  + 

+ 

0.60  + 

0.50  + 

0.40+ 

V 

V 

+ 

- - 

V  — - 

-v 

+ 

0.30  -t" - 

V 

+ 

0.20  + 

+ 

0  .  10  V 

+ 

0.00  +....+. 

...+....+... 

.  +  . 

.  .  .  +  . 

...  +  .. 

.  .+  . 

....  +  ....+  .  .  . 

+ 

0  .  1 

.2  .3 

.4 

.5 

.6 

.7 

.8  .9  1  . 

0 

THE  ABSCISSA 

IS  Z 

MIN  = 

1  .  OOE+00 

MAX=  2.31E+01 

Figure  C.27.  V  vs  Z-Axis  At  IZ=21  To  IZ=26 
For  IX=I  And  IY=13 


PATCH(PL0T16  ,PROFIL,  1,  5,  13,  13, 

PL0T(PL0T16  ,H1  .  5.000E+05,  6.000E+05) 

TABULATION  OF  ABSCISSA  AND  ORDINATES... 


1,  1) 


2.500E-01 
7 .500E-01 
1 .250E+00 
1 .750E+00 
2.250E+00 
VARIABLE 
MINVAL= 
MAXVAL= 
C£LLAV= 


l.OO  +. 
0.90  + 
0.80  + 
0.70  H_ 
0.60  + 
0,50  + 
0.40  + 
0.30  + 
0.20  + 
0.10  + 
0.00  +. 
0 

THE  ABt 


HI 

5.653E+05 

5.647E+05 

5.638E+05 

5.628E+05 

5.620E+05 

HI 

5.000E+05 
6  .  OOOE+05 
5.637E+05 


.1  .2  .3  .4  .5  .6  .7  .8  .9  1.0 

5CISSA  IS  X  .  MIN=  2.50E-01  MAX=  2.25E+00 


Figure  C.28.  H  vs  X-Axis  At  IX=1  To  1X=5 
For  IY=13  And  IZ=23 


115 


PATCH(PL0T17  ,PROFIL,  1,  5,  13,  13,  23,  23, 

PL0T(PL0T17  ,V1  ,  l.OOOE+01,  4.500E+01) 


TABULATION 

OF  ABSCISSA  AND  ORDINATES 

X 

VI 

2.500E-01 

2.223E+01 

7.S00E-01 

2.306E+01 

1 .250E+00 

2.448E+01 

1 .75QE+00 

2.591E+01 

2.250E+00 

2.6'ilE-hOl 

VARIABLE 

VI 

MINVAL= 

1 . OOOE+01 

MAXVAL= 

9.500E+01 

CELLAV= 

2.9A2E+01 

1) 


1.00  _ _ 

0,90  +  + 
0.80+  + 
0.70  +  + 
0.60  +  _ ^ _  + 

olso  V  + 
0.20  +  + 
0,10  +  + 
0.00 

0  .1  .2  .3  .9  .5  .6  .7  .8  .9  1.0 
THE  ABSCISSA  IS  X  .  MIN=  2.50E-01  MAX=  2.25E+Q0 


Figure  C.29.  V  vs  X-Axis  At  IX=1  to  IX=5 
For  IY=i3  And  IZ=23 


lie 
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